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>< ABSTRACT 

S3 

We present high angular resolution 1.3 mm and 850 /im dust continuum data 
obtained with the Submillimeter Array toward 33 Class protostars in nearby 
clouds (distance < 500 pc), which represents so far the largest survey toward pro- 
tostellar binary/multiple systems. The median angular resolution in the survey 
is 2'.'5, while the median linear resolution is approximately 600 AU. Compact 
dust continuum emission is observed from all sources in the sample. Twenty-one 
sources in the sample show signatures of binarity/multiplicity, with separations 
ranging from 50 AU to 5000 AU. The numbers of singles, binaries, triples, and 



quadruples in the sample are 12, 14, 5, and 2, respectively. The derived multi- 
plicity frequency (MF) and companion star fraction (CSF) for Class protostars 
are 0.64 ±0.08 and 0.91 ±0.05, respectively, with no correction for completeness. 
The derived MF and CSF in this survey are approximately two times higher than 
the values found in the binary surveys toward Class I young stellar objects, and 
approximately three (for MF) and four (for CSF) times larger than the values 
found among main sequence stars, with a similar range of separations. Further- 
more, the observed fraction of high order multiple systems to binary systems in 
Class protostars (0.50 ±0.09) is also larger than the fractions found in Class! 
young stellar objects (0.31 ±0.07) and main sequence stars (<0.2). These re- 
sults suggest that binary properties evolve as protostars evolve, as predicted by 
numerical simulations. The distribution of separations for Class protostellar 
binary/multiple systems shows a general trend in which companion star frac- 
tion increases with decreasing companion separation. We find that 67% ± 8% of 
the protobinary systems have circumstellar mass ratios below 0.5, implying that 
unequal-mass systems are preferred in the process of binary star formation. We 
suggest an empirical sequential fragmentation picture for binary star formation, 
based on this work and existing lower resolution single-dish observations. 

Subject headings: binaries: general — ISM: clouds — ISM: dust, extinction — 
stars: formation — techniques: interferometric 



1. INTRODUCTION 

Over the past two decades, our knowledge of the formation and evolution of low-mass 
stars has made significant progress (see, e.g., Reipurth et al. 2007 for several reviews). It is 
widely accepted that low-mass stars form from the gravitational collapse of molecular cloud 
cores (see, e.g., Shu et al. 1987; McKee & Ostriker 2007). Initially, these cores, generally 
referred to as prestellar cores, are cold dense condensations with infall motions that have 
not yet formed central stellar objects (Andre et al. 2000; 2009). Resulting from the collapse 
of prestellar cores. Class objects are the youngest accreting protostars observed right after 
point mass formation, when most of the mass of the system is still in the surrounding 
dense core/envelope (Andre et al. 2000). As the protostellar envelope dissipates through 
accretion and ejection of circumstellar material. Class protostars evolve into Class I young 
stellar objects (YSOs; Lada 1987). Class I YSOs subsequently evolve into pre-main sequence 
(PMS) stars surrounded by circumstellar disks, but without dense envelopes. Eventually, 
PMS stars evolve into main sequence (MS) stars. 



One of the major puzzles in our understanding of star formation concerns the origin 
of binary stars (Tohhne 2002). A series of surveys toward solar- type MS stars have found 
that about 40%-60% of MS stars are actually binary systems, with periods from less than 
one day to over 10 million years (e.g., Duquennoy & Mayor 1991; Fischer & Marcy 1992). 
The binary frequency found among PMS stars is even higher than that among the MS stars 
(e.g., Mathieu 1994). For example, in the survey toward PMS stars in the Taurus star 
formation region, Kohler & Leinert (1998) found that for systems with separations from 
18 AU to 1800 AU the binary frequency is 48.9% ±5.3%, which is larger than a factor of 
two compared to the binary frequency in solar-type MS stars with the same separation 
range. Recent surveys toward Class! YSOs show that binary systems are also common in 
the Class! phase (see, e.g., !!aisch et al. 2004; Duchene et al. 2004, 2007; Connelley et 
al. 2008a, b). For example, the large survey by Connelley et al. (2008a, b) found that for 
Class! YSOs the binary frequency for separations from 100 AU to 4500 AU is about 43%. 
The high frequency of binary systems found in different stellar evolution stages has led to 
a widely held opinion that most stars form in binary/multiple systems (e.g., Mathieu 1994; 
Goodwin et al. 2007). !n contrast, however, Lada (2006) argued that the most common 
outcome of the star formation process is a single star rather than a binary star, based on two 
observational facts: (1) that most low-mass stars (stellar mass <O.5M0) are single objects 
(e.g., Delfosse et al. 2004; Siegler et al. 2005), and (2) that these low-mass stars represent 
the majority of the field stars (Kroupa 2002; Chabrier 2005). 

Do most stars form as singles or binaries? Before trying to answer this question still in 
debate, it must be noted that multiple stellar systems undergo dynamical evolution at the 
beginning of their formation phase (e.g., Reipurth & Clarke 2001; Reipurth et al. 2010). Dy- 
namical evolution will modify the binary population by reducing the overall binary fraction 
and by hardening the remaining binaries on a timescale of less than 10^ yr (see Goodwin et 
al. 2007). Therefore, neither observations of PMS/Class! young stars nor observations of 
field dwarfs can constrain the original properties of binary/multiple systems. Clearly, the 
most effective way to answer this question is to extensively observe the earliest stages of the 
star formation process. 

Class objects, which represent the youngest protostars at an age of a few 10^—10^ yr 
(Andre et al. 2000; Evans et al. 2009), provide an important opportunity to probe the earliest 
phase of the star formation process, where most of the initial information of (binary) star 
formation is expected to still be preserved. To understand the formation of binary stars, 
it is therefore critical to observe large samples of Class protostars in nearby molecular 
clouds to answer a number of key questions, e.g., how common is binarity/multiplicity in 
the Class protostellar phase? What are the initial distributions of separations and mass 
ratios? Unfortunately, direct observations of the Class protostellar stages, when the main 



collapse has started but no optical or near-infrared emission emerges from the protostar 
through the opaque infalling envelope, were long hampered by the low angular resolution 
of single-dish (sub-)millimeter telescopes. Only the high angular resolution attained with 
current millimeter interferometers enables us to directly observe the formation phase of 
binary stars. In the past decade, a handful of Class protostellar binary (protobinary) 
systems have been found using different interferometric arrays (e.g., Looney et al. 2000; 
Reipurth et al. 2002; Launhardt 2004), and there are increasing interferometric studies of 
binarity in the protostellar phase (e.g., Volgenau et al. 2006; Maury et al. 2010). However, 
these studies were limited to small samples of objects and no statistically strong conclusions 
could be drawn from them. 

On the theoretical side, numerous simulations support the hypothesis that the fragmen- 
tation of molecular cloud cores is the main mechanism for the formation of binary/multiple 
stellar systems (see reviews by Bodenheimer et al. 2000, Tohline 2002, and Goodwin et 
al. 2007). Yet, many key questions concerning this fragmentation process (e.g., exactly 
when, where, why, and how) are still under debate. In general, fragmentation models can 
be classified depending on when the fragmentation takes place during the different collapse 
phases of the coraM Initial clump/filament fragmentation takes place prior to the collapse of 
individual cores, and results in cores with separations of ~ 10^— 10"^ AU (see, e.g., Kauffmann 
et al. 2008; Launhardt et al. 2010), while prompt (isothermal) fragmentation takes place 
at the end of the core's isothermal collapse phase, and results in fragments with separations 
of ~10^— lO'^AU (see Tohline 2002). On the other hand, adiabatic fragmentation occurs 
when the collapse evolves into the adiabatic phase, and results in fragments with separations 
of ~3— 300AU (Machida et al. 2005; 2008). Finally, secondary fragmentation takes place 
during the core's secondary collapse phase, and results in fragments with separations of 0.01 
to 0.1 AU (Machida et al. 2008). In all these cases, fragmentation takes place when there is 
some angular momentum in the core (see the review by Goodwin et al. 2007). The source of 
the angular momentum may be ordered rotational motions (i.e., rotational fragmentation) 
or turbulence (i.e., turbulent fragmentation). 

These different possible models of fragmentation result in different binary/multiple sys- 
tems with a wide range of properties. These properties are largely determined by the accre- 
tion process which, in turn, strongly depends on the initial conditions of the cloud cores, e.g., 
the initial distributions of mass and angular momentum (see Bate & Bonnell 1997). It must 
be noted that the constraints used on current binary formation models mostly come from the 
observations of MS and PMS stars. These observations cannot provide information on the 



"'^See, e.g., Figure 2 in Tohline 2002 or Figure 9 in Andre et al. 2009 for a discussion of the different 
evolutionary phases of core collapse. 



initial conditions of fragmenting cores. Hence, it is critical to observe the early phase of core 
fragmentation and study in detail their properties (e.g., mass and separation distributions), 
in order to put direct constraints on different fragmentation models. 

To achieve a comprehensive understanding of binary star formation, we have started a 
systematic program to observe at different (sub-) millimeter interferometers a large sample 
of low-mass prestellar and protostellar cores. The early survey was conducted at the Owens 
Valley Radio Observatory (OVRO) millimeter array (Chen et al. 2007; R. Launhardt et 
al. in preparation), and was then continued with the Australia Telescope Compact Array 
(ATCA; Chen et al. 2008a), the IRAM Plateau de Bure Interferometer (PdBI; Chen et al. 
2009; X. Chen et al. in preparation), and the Submillimeter Array (SMA; Chen et al. 2008b; 
Chen & Arce 2010). As part of this program, we present in this paper high angular resolu- 
tion (sub-) millimeter continuum observations toward 33 Class protostars using the SMA, 
which represents so far the largest survey of protobinary systems. High angular resolution 
(sub-) millimeter continuum emission traces the circumstellar dust from the inner envelopes 
and disks surrounding protostars, which we used to resolve binary /multiple systems, and 
measure their masses, sizes, and spatial distribution. In Section 2 we describe the target 
list, observations, and data reduction. Observational results are presented in Section 3 and 
discussed in Section 4. The main conclusions of this study are summarized in Section 5. We 
give a detailed description of the individual sources in Appendix A. 



2. SAMPLE, OBSERVATIONS, AND DATA REDUCTION 

2.1. Sample 

The sources in this study were collected from the SMA science archive. A total of 33 
sources were selected by using the following three criteria: (1) available in the SMA archivaj 
(2) confirmed classification as a Class object by other studies; and (3) known distance less 
than 500 pc. These sources were observed at the SMA between 2004 and 2009 with different 
science goals (e.g., for studying outflow, kinematics, or chemistry in the Class protostars). 
Only for a few sources (i.e., IRAS 03282, NGC 1333 SVS 13, and CG 30), the proposed science 
goal was to study the formation of binary stars. The target list and their basic properties 
(i.e., bolometric luminosity Lboi, bolometric temperature Tboi, and envelope mass Menv) 
are summarized in Table 1. Most sources in Table 1 were cataloged by both Andre et al. 

^All the SMA data used in this work are pubhcly available from the SMA archive. Some of the data 
are unpublished and in Tables 2-3 we list the name of the principle investigator (PI) responsible for the 
observations. For data that had been published we list the paper that first presented the data in Tables 2-3. 



(2000) and Froebrich (2005). We refer readers to these two papers and references therein 
for more details about the basic properties of these sources. Note that several sources may 
represent transition objects between the Class and Class! phases (see Froebrich 2005 for 
more detailed discussion). 

Distances to the different sources were obtained from the literature. Early estimates of 
the distance to Perseus indicated that this region is 350 pc away (e.g., Herbig & Jones 1983). 
However, more recent studies indicate that this region is significantly closer. Photometric 
distances towards the western part of Perseus place this cloud at 220 pc (Cernis 1990; Cernis 
& Straizys 2003), consistent with recent VLBI observations that give a distance to NGC 1333 
of 240±20pc (Hirota et al. 2008, 2011). We therefore adopt a distance of 240 pc for all 
objects in Perseus, in agreement with other recent studies of this cloud (e.g., Enoch et al. 
2006; Evans et al. 2009; Arce et al. 2010, 2011). The basic properties (bolometric luminosity 
and envelope mass) of the Perseus sources listed in Andre et al. (2000) and Froebrich (2005) 
were modified using the new distance of 240 pc (see Table 1). For sources in the Taurus star 
forming region, we adopt a distance of 140 pc (Elias 1978; Kenyon et al. 1994; Loinard et 
al. 2007b). For sources in the Orion region, we adopt a distance of 440 pc, obtained by the 
VLBI observations of Hirota et al. (2007). For sources in the Ophuichus cloud, we use a 
distance of 120 pc, derived from recent VLBA observations by Loinard et al. (2008). For 
the Lupus region, the distance is estimated to be 130 pc (Murphy et al. 1986). For isolated 
objects, the distance references are as follows: CB 17 (Launhardt et al. 2010), CG30 (Knude 
et al. 1999), L483 (Dame & Thaddeus 1985), L723 (Goldsmith et al. 1984), B335 (Tomita 
et al. 1979), L1157 (Straizys et al. 1992), and L1251B (Kun & Prusti 1993). 



2.2. Observations 



Observations were carried out with the SMAJj (Ho et al. 2004) between 2004 and 2009. 
Five different array configurations, Subcompact (Subcom), Compact (Com), Compact-North 
(Com-N; for several equatorial and southern sources). Extended (Ext), and Very Extended 
(V-Ext), were used in the observations, with baselines ranging from about 10 m to 500 m. 
Only one source, IRAS 16293-2422, was also observed with the Extended Submillimeter 
Array (eSMA^). 



■'The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the 
Academia Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and 
the Academia Sinica. 

*See 



http://www.jach.hawaii.edu/JCMT/eSMA/ 



for more details. 



Table 1. Basic properties of target sources 



Source 


Other (or short) 


Association 


Distance 


Lhoi 


Tboi 


Menv 


Class 


References'' 


Name 


Name" 


Region 


[pc] 


[L0] 


[K] 


[Mq] 






L1448 IRS2 


IRAS 03222+3034 


Perseus 


240 


1.7-4.1 


43-63 


0.6-1.4 





1, 2, 3, 4 


L1448 N 


IRAS 03225+3034 


Perseus 


240 


4.3-17 


53-90 


1.5-4.5 


O/I 


1, 2, 3, 4 


L1448 C 


L1448-mm 


Perseus 


240 


4.4-6.3 


49-69 


0.7-1.9 





1, 2, 3, 4 


PER 065 


Per-Bolo 30 


Perseus 


240 


0.03-0.22 


31-37 


0.4-2.1 





3, 4 


NGC 1333 IRAS 2A 


IRAS 2A 


Perseus 


240 


19-23 


50-69 


1.0-2.8 





1, 2, 3, 4 


NGC1333IRAS2B 


IRAS 2B 


Perseus 


240 


19-23 


50-69 


1.0-2.8 





1, 2, 3, 4 


NGC 1333 SVS 13 


SVS 13 


Perseus 


240 


7-33 


30-180 


1.3-3.3 


O/I 


1, 2, 3, 4 


NGC 1333 IRAS 4A 


IRAS 4A 


Perseus 


240 


4.2-8.5 


34-51 


2.7-7.8 





1, 2, 3, 4 


NGC 1333 IRAS 4B 


IRAS 4B 


Perseus 


240 


1.6-7.6 


36-54 


1.4-3.7 





1, 2, 3, 4 


IRAS 03282+3035 


IRAS 03282 


Perseus 


240 


0.7-1.2 


33-60 


0.4-1.0 





1, 2, 3, 4 


Per Bl-c 


Bl-c 


Perseus 


240 


1.8-3.7 


53-76 


2.1-18 


O/I 


3, 4, 5, 6 


Per Bl-b 


Bl-b 


Perseus 


240 


1.2-2.5 


<25 


2.8-26 





3, 5, 7 


HH211-MM 




Perseus 


240 


1.5-4.5 


24-31 


0.5-3.0 





1, 2, 3, 4 


IC 348 MMS 




Perseus 


240 


0.5-1.9 


49-59 


1.8-10 





3, 4 


CB 17 IRS 


L1389 


Isolated 


250 


0.6 


>55 


<4 


O/I 


8 


IRAM 04191+1522 


IRAM 04191 


Taurus 


140 


0.12-0.28 


18-36 


0.5 





1, 2, 9, 10 


L1521F-IRS 


L1521F/MC27 


Taurus 


140 


0.1-0.36 


20 


0.7-4.8 





11, 12, 13 


L1527 IRS 


IRAS 04368+2557 


Taurus 


140 


1.9-2.2 


36-60 


0.4-0.9 


O/I 


1, 2, 14 


HH 114 MMS 




Orion 


440 


<25 


40-84 


1.4-2.8 





1, 2 


OMC3 MM 6 




Orion 


440 


<60 


30-72 


6-12 





1, 2 


NGC2024-FIR5 




Orion 


440 


>10 


20 


8.2-15 





1, 2 


HH212-MM 




Orion 


440 


7.7-14 


<56 


0.3-1.2 





1, 2 


HH 25 MMS 




Orion 


440 


6.0-7.2 


34-41 


0.5-1.2 





1, 2 


CG30 


BHR12 


Isolated 


200 


10-17 


74-117 


1.1-2.3 


O/I 


1, 2 


B228 


IRAS 15398-3359 


Lupus 


130 


0.9-1.2 


48-61 


0.3-0.4 





2, 14 


VLA 1623 




Ophuichus 


120 


0.4-1.2 


<59 


0.2-0.5 





1, 2, 4 


Oph MMS 126 


IRAS 16253-2429 


Ophuichus 


120 


0.25 


35 


0.51 





4 


IRAS 16293-2422 


IRAS 16293 


Ophuichus 


120 


6.9-14 


41-54 


1.2-2.8 





1, 2, 4 


L483 


IRAS 18148-0440 


Isolated 


200 


9-13 


50-54 


0.3-1.8 


O/I 


1, 2, 14 


L723VLA2 


IRAS 19156+1906 


Isolated 


300 


3-3.3 


47-50 


0.6-1.6 





1, 14 


B335 


IRAS 19345+0727 


Isolated 


250 


3.0-3.1 


28-45 


0.7-1.2 





1, 2, 14 


L1157 


IRAS 20386+6751 


Isolated 


325 


5.8-11 


42-60 


0.5-2.6 





1, 2, 14 


L1251B IRSl 


IRAS 22376+7455 


Isolated 


300 


<10 


83-87 


2.0-2.1 


O/I 


2, 15, 16 



''The IRAS source is not always associated with the mm source. 

''References for the basic properties of the objects. — (1) Andre et al. (2000); (2) Frobcrich (2005); (3) Hatchell ct al. (2007); 
(4) Enoch et al. (2009); (5) Matthews & Wilson (2002); (6) Matthews et al. (2006); (7) Hirano et al. (1999); (8) Launhardt 
et al. (2010); (9) Andre et al. (1999); (10) Dunham et al. (2006); (11) Crapsi et al. (2004); (12) Bourke et al. (2006); (13) 
Terebey et al. (2009); (14) Shirley et al. (2000); (15) J. Lee et al. (2006); (16) J. Lee et al. (2007). 
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The 230 GHz and 345 GHz receivers were used in the observations. Both receivers have 
two sidebands (lower and upper), with 2 GHz bandwidth each (separated by 10 GHz). The 
receiver setups were typically tuned to the CO (2-1) line in the 230 GHz observations (see, 
e.g., Chen et al. 2008b; 2010) and to the CO (3-2) line in the 345 GHz observations (see, 
e.g., J0rgenson et al. 2007; Lee et al. 2009). The total continuum bandwidth, combining the 
line-free portions of the two sidebands, are about 3.0-3.8 GHz for the 1.3 mm and 850 /xm 
dust continuum. We present only the continuum results in this paper. At 230 GHz, system 
temperatures generally range from 80 to 300 K (depending on elevation), with a typical value 
of 150 K. At 350 GHz, system temperatures range from 150 to 500 K, with a typical value of 
250 K. The SMA primary beams are ~ 55" at 230 GHz and ~ 36" at 345 GHz. 

The visibility data were reduced using the MIR package (Qi 2005). During the reduction, 
visibility points with clearly deviating phases and/or amplitudes were flagged. The passband 
(spectral response) was calibrated through observations of available planets (e.g., Saturn) 
and strong quasars (e.g., 3C 454.3) at the beginning and/or end of each track. Time variation 
of amplitude and phase was calibrated through frequent observations of quasars near to each 
source. The flux density scale was calibrated by observing planets such as Uranus and Ceres. 
For a few tracks, strong quasars (e.g., 3C273) with stable flux densities during the observing 
dates, were also used as flux calibrators. We estimate a typical flux accuracy of about 20%- 
30%. The calibrated visibility data were then imaged using the MIRIAD toolbox (Sault 
et al. 1995). Each continuum map was deconvolved down to the la theoretical rms noise 
level using the MIRIAD clean routine. The synthesized beam sizes (with robust weighting 
1.0) and rms noise levels are summarized in Tables 2 and 3. Figure 1 shows the beam sizes 
plotted as a function of the corresponding linear resolutions. The median angular resolution 
in this survey is 3"1 at A 1.3 mm and 1"0 at A850/im, corresponding to the median linear 
resolutions of approximately 640 AU and 220 AU, respectively. Combining the observations 
at the two wavelengths, the median angular and linear resolutions in this survey are 2"5 and 
600 AU, respectively. Figure 2 shows the distribution of the rms noise levels. The median 1 a 
theoretical noise is l.OmJybeam"^ at Al. 3 mm and 1.7mJybeam^^ at A850/im, while in the 
cleaned images, the median la noise is 2.3mJybeam~^ and 5.5mJybeam~^, respectively. 



3. RESULTS: SMA SUBMILLIMETER AND MILLIMETER MAPS 

Figures 3-17 present the SMA submillimeter and millimeter maps for the sources listed 
in Table 1. Most sources in Table 1 were observed with SCUBA on the James Clerk Maxwell 
Telescope (JCMT). Complementary SCUBA submillimeter images are also shown for selected 
regions and sources, to illustrate the relative locations of the protostellar cores within the 
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clouds and/or the overall morphology of the sources. The SCUBA data for the sources in 
the Perseus and Ophuichus clouds are available from the COMPLETE databasqj The data 
for other objects were retrieved from the JCMT science archivqj 

In the SMA observations, millimeter and/or submillimeter dust continuum emission is 
detected from all the sources in the sample. It must be noted that in this work we regard a 
local density peak or enhancement as a source if its emission is > 5 a levels with respect to 
its immediate surrounding. Except several very low luminosity sources such as PER 065 (see 
Figure 4), L1521F-IRS (see Figure 9c), and Oph-MMS 126 (see Figure 15c), most sources in 
the sample show centrally peaked emission in the SMA dust continuum images, suggesting 
the existence of a compact inner envelope and/or circumstellar disk in these sources. 

As shown in the SCUBA images (angular resolution of ~ 14" at 850 /zm), most sources 
in the sample appear to be single cores (except CG30). At the high angular resolution 
obtained with the SMA (median angular resolution of 3'.'1 at 1.3 mm and of I'.'O at 850 /im), 
22 sources in the sample show signatures of binarity/multiplicity, with angular separations 
from 0'.'3 (HH211 MMS) to 22" (CG30). Detailed descriptions of the individual sources are 
presented in Appendix A. 

Tables 4 and 5 list the derived parameters (e.g., positions, integrated fluxes, and decon- 
volved FWHM sizes) for each of the sources/components (> 5 a level) from circular Gaussian 
fits to the visibilities. The statistical uncertainties on the derived parameters are generally 
small and are mainly determined by the rms noise levels of the images (see Tables 2 and 
3). On the other hand, the calibration uncertainties of the fluxes are 20%-30%, and thus 
dominate the uncertainties of the derived fluxes. For a few faint components in the systems 
such as IRAS 03282, NGC 1333 IRAS2A, and HH212 MMS, their fluxes are uncertain to 
within a factor of two, due to uncertainties in separating their emission from the surrounding 
envelope emission. 

Figure 18 shows the distribution of fluxes derived in this work. The 1.3 mm fluxes 
range from 1.1 mJy to 3.0 Jy, with a median value of 100 mJy, while the 850 /xm fluxes range 
from 20mJy to 3.2 Jy, with a median value of 430 mJy. Compared with the total fluxes 
obtained from the single-dish IRAM-30m images at 1.3 mm (see, e.g., Kauffmann et al. 
2008; Launhardt et al. 2010) and the SCUBA images at 850 /im (see, e.g.. Chandler & Richer 
2000; Shirley et al. 2000; Sandell &; Knee 2001), the SMA resolved out a significant fraction 



^See Ridge et al. (2006) and http://www.cfa.harvard.edu/COMPLETE/ 



for more details. 



^http://cadcwww.dao.nrc.ca/jcnit/. The JCMT archive at the Canadian Astronomy Data Center is 
operated by the Herzberg Institute of Astrophysics, National Research Council of Canada. 
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(generally 80%-90%) of the extended emission and detected only the compact emission in 
the center of protostellar cores (i.e., inner envelopes and circumstellar disks). 

Assuming that the 1.3 mm and 850 /im continuum emission is optically thin, the total 
gas mass (Mgas) of the individual sources can be derived from the (sub-) millimeter fluxes 
using the formula: 

F D"^ 

where F^ is the flux of the source, D is the distance to the source, Tj is the dust temperature, 
/tm(^) is the dust opacity (mass absorption coefficient per gram of dust, e.g., 1.75 cm^g^"*^ at 
850 /im for 0H5 dust; see Ossenkopf & Henning 1994), and By{T^ is the Planck function at 
frequency v and dust temperature T^. In practice, we use the same method as in J0rgensen 
et al. (2007) to calculate the masses (adopting a dust temperature of 30 K for all sources), 

Mi.3^^ = 1.3 Mo m^\ f;^V X {exp[0.36 (^^\ - 1}, (2) 



IJy y V200 pc 



0.85 mill 



D V . , /30 K 



Mo.85mm = 0.18 Me [-^yW) V200^J "" {^M^-^^ \^ T " ^^- ^^^ 

Column 5 of Tables 4 and 5 lists the total gas masses for each sources, calculated from the 
1.3 mm and 850 /im fluxes. The relative uncertainty of the derived masses due to calibration 
errors of the (sub-) millimeter fluxes is about 20%-30% for all sources. Moreover, due to 
the uncertainties in the dust opacity and temperature, the masses obtained this way may 
contain uncertainties of up to a factor of two. 

Figure 18 shows the distribution of gas masses for the objects in this sample. The masses 
derived from the 1.3 mm observations range from 0.001 M0 to 2AMq, with a median value 
of 0.06 Mq, while the masses derived from the 850 /xm observations range from 0.005 M© to 
1.7 Mq, with a median value of 0.12 Mq. Compared with the masses derived from single-dish 
observations, e.g., a median value of ~ 0.5 Mq derived from IRAM-30m 1.3 mm observations 
for protostars (see Kauffmann et al. 2008), the masses derived from this SMA survey are 
relatively small and represent only lower limits to the core masses. As indicated above, 
most of the (sub-) millimeter continuum emission is resolved out in the SMA images, and it 
is very likely that detected emission arises from the inner envelope and circumstellar disk. 
Therefore, the derived masses may also be regarded as upper limits to the disk masses (see, 
e.g., J0rgensen et al. 2009 for further discussion). 
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Table 4. SMA 1.3 mm dust continuum results. 



Source 


R.A.° 


Dec." 


5^" 


FWHM" 


Mgas'' 




(J2000) 


(J2000) 


[mjy] 


(maj." X mill.") 


[Mq] 


L1448 IRS2 


03:25:22.38 


+30:45:13.3 


180+36 


3.9x2.3 


0.15 


L1448NA 


03:25:36.48 


+30:45:21.8 


140+30 


3.0x1.4 


0.11 


L1448NB 


03:25:36.33 


+30:45:14.8 


660+130 


2.3x1.6 


0.54 


L1448NNW 


03:25:35.66 


+30:45:34.2 


70+14 


3.5x2.5 


0.057 


L1448CN 


03:25:38.88 


+30:44:05.5 


200+40 


1.6x0.9 


0.16 


L1448CS 


03:25:39.05 


+30:43:59.3 


20+5 


5.7x1.9 


0.016 


PER 065 A 


03:28:39.17 


+31:05:55.6 


4.4+1 


8.1x6.5 


0.004 


PER 065 B 


03:28:38.86 


+31:05:56.8 


1.1+0.2 


2.3x1.3 


0.001 


PER 065 C 


03:28:38.96 


+31:06:04.2 


1.5+0.3 


2.3x2.1 


0.001 


NGC 1333 IRAS 2A 


03:28:55.58 


+31:14:37.2 


330+66 


1.8x1.3 


0.27 


NGC 1333 IRAS 2B 


03:28:57.38 


+31:14:15.8 


50+10 


3.8x2.0 


0.04 


NGC 1333 SVS 13 A 


03:29:03.75 


+31:16:03.6 


450+90 


2.7x2.0 


0.37 


NGC1333SVS13VLA3 


03:29:03.41 


+31:16:01.5 


22+5 


(unresolved) 


0.02 


NGC 1333 SVS 13 B 


03:29:03.07 


+31:15:51.8 


460+90 


2.8x2.3 


0.37 


NGC 1333 SVS 13 C 


03:29:01.97 


+31:15:38.1 


50+10 


2.2x0.7 


0.04 


NGC 1333 IRAS 4A 


03:29:10.51 


+31:13:31.3 


3000+600 


2.4x1.7 


2.43 


NGC 1333 IRAS 4B 


03:29:12.01 


+31:13:08.2 


920+180 


1.4x0.9 


0.75 


NGC 1333 IRAS 4C 


03:29:12.84 


+31:13:07.0 


250+50 


0.7x0.6 


0.20 


IRAS 03282 


03:31:20.94 


+30:45:30.2 


290+60 


1.4x1.3 


0.24 


Pcr-Bl-c 


03:33:17.87 


+31:09:32.3 


300+60 


3.1x2.4 


0.24 


Per-Bl-bN 


03:33:21.19 


+31:07:43.6 


240+50 


4.6x3.2 


0.19 


Per-Bl-bS 


03:33:21.33 


+31:07:26.3 


340+70 


3.2x1.8 


0.28 


IC348MMS1 


03:43:57.05 


+32:03:05.0 


260+50 


2.5x1.9 


0.21 


IC348 MMS2 


03:43:57.73 


+32:03:10.5 


65+15 


4.4x2.6 


0.053 


CB17IRS 


04:04:33.76 


+56:56:16.5 


6.3+1.3 


3.4x2.5 


0.006 


CB17MMS 


04:04:35.78 


+56:56:03.4 


3.6+0.7 


3.6x1.1 


0.003 


IRAM 04191 IRS 


04:21:56.89 


+15:29:46.1 


4.5+1.5 


2.8x2.7 


0.001 


IRAM 04191 MMS 


04:21:56.36 


+15:29:48.9 


2.5+1.0 


3.7x1.1 


< 0.001 


L1521F IRS 


04:28:38.91 


+26:51:34.6 


15+5 


5.0x2.1 


0.004 


L1527 


04:39:53.90 


+26:03:09.8 


190+40 


1.0x0.5 


0.052 


HH 114 MMS 


05:18:15.24 


+07:12:02.7 


490+100 


1.9x1.6 


1.34 


HH 212 MMS 


05:43:51.41 


-01:02:53.1 


76+15 


0.6x0.5 


0.21 


CG30N 


08:09:33.13 


-36:04:58.3 


260+50 


3.3x1.9 


0.15 


CG30S 


08:09:32.73 


-36:05:19.3 


100+20 


4.6x2.4 


0.06 


B228 


15:43:02.20 


-34:09:07.6 


150+30 


5.5x3.9 


0.036 


B228(l.lmm) 


15:43:02.13 


-34:09:07.3 


135+30 


5.4x4.3 


0.032 


VLA 1623 W 


16:26:25.64 


-24:24:29.3 


42+10 


1.1x0.7 


0.009 


VLA1623EA 


16:26:26.39 


-24:24:30.7 


79+15 


1.0x0.8 


0.016 


VLA 1623 EB 


16:26:26.31 


-24:24:30.4 


71+15 


0.6x0.6 


0.014 


Oph-MMS 126 


16:28:21.58 


-24:36:23.6 


11+2 


3.6x1.7 


0.002 


IRAS 16293 A 


16:32:22.87 


-24:28:36.6 


620+120 


1.7x1.4 


0.13 


IRAS 16293 B 


16:32:22.62 


-24:28:32.6 


1300+260 


1.6x1.4 


0.26 


L723VLA2MMS1 


19:17:53.68 


+19:12:19.8 


59+12 


2.7x2.0 


0.075 


L723VLA2MMS2 


19:17:53.91 


+19:12:18.4 


58+12 


3.2x1.4 


0.074 


B335 


19:37:00.93 


+07:34:09.9 


145+30 


4.5x2.4 


0.13 
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Table 4 — Continued 



Source 


R.A.'^ 


Dec.'' 


S.'' 


FWHM" 


Mga,6 




(J2000) 


{J2000) 


[mjy] 


(maj." X min.") 


[Mq] 


L1157 


20:39:06.18 


+68:02:16.3 


280+56 


3.1x3.0 


0.42 


L1251BIRS1 


22:38:46.68 


+75:11:33.1 


48+10 


4.2x2.3 


0.061 


L1251BIRS2 


22:38:52.94 


+75:11:23.2 


21+4 


3.6x3.3 


0.027 


L1251BSMAN 


22:38:48.95 


+75:11:33.4 


17+4 


4.9x4.7 


0.022 


L1251BSMAS 


22:38:49.09 


+75:11:25.0 


28+6 


6.1x4.2 


0.035 



''Center position, flux, and FWHM size of the dust continuum sources derived from 
Gaussian fitting. 

''Total gas mass derived from the SMA 1.3 mm dust continuum observations. Note 
that the gas mass is from compact emission only, not the whole envelope. 



Table 5. SMA 850 fira dust continuum results. 



Source 


R.A." 


Dec.'' 


S.^ 


FWHM" 


Mga," 




(J2000) 


{J2000) 


[mJy] 


(maj." X min.") 


[Mo] 


L1448CN 


03:25:38.87 


+30:44:05.3 


410+80 


0.6x0.4 


0.078 


L1448CS 


03:25:39.14 


+30:43:60.0 


40+10 


1.2x0.6 


0.008 


NGC 1333 IRAS 2A1 


03:28:55.57 


+31:14:37.1 


800+160 


1.1x0.8 


0.15 


NGC 1333 IRAS 2A2 


03:28:55.55 


+31:14:38.5 


40+10 


(unresolved) 


0.008 


NGC 1333 IRAS 4A1 


03:29:10.53 


+31:13:31.0 


3200+600 


1.4x0.4 


0.61 


NGC 1333 IRAS 4A2 


03:29:10.44 


+31:13:32.2 


1600+300 


1.1x0.1 


0.30 


IRAS 03282 MMSl 


03:31:20.94 


+30:45:30.3 


340+70 


0.2x0.1 


0.065 


IRAS 03282 MMS2 


03:31:21.08 


+30:45:29.9 


30+10 


(unresolved) 


0.006 


HH211 SMMl 


03:43:56.80 


+32:00:50.3 


80+20 


0.1x0.1 


0.015 


HH211 SMM2 


03:43:56.79 


+32:00:50.0 


25+10 


<0.1 


0.005 


HH114MMS 


05:18:15.23 


+07:12:02.5 


1100+220 


1.3x1.2 


0.70 


OMC3MMS6 SMMl 


05:35:23.43 


-05:01:30.6 


1350+200 


1.0x0.8 


0.86 


OMC3MMS6 SMM2 


05:35:23.52 


-05:01:41.2 


550+100 


1.8x1.7 


0.35 


NGC2024 FIR5 SMMl 


05:41:44.26 


-01:55:40.9 


2640+500 


2.9x1.8 


1.69 


NGC2024 FIR5 SMM2 


05:41:44.51 


-01:55:42.4 


1030+200 


3.5x3.0 


0.66 


HH 212 MMSl 


05:43:51.41 


-01:02:53.1 


260+50 


2.0x1.0 


0.17 


HH212MMS2 


05:43:51.38 


-01:02:53.0 


20+4 


(unresolved) 


0.013 


HH212MMS3 


05:43:51.48 


-01:02:53.8 


27+5 


0.5x0.4 


0.018 


HH 25 SMMl 


05:46:07.26 


-00:13:30.7 


510+100 


2.8x2.4 


0.33 


HH25SMM2 


05:46:07.34 


-00:13:43.4 


430+80 


1.3x0.7 


0.27 


HH25SMM3 


05:46:07.58 


-00:13:53.7 


150+30 


4.8x4.6 


0.10 


IRAS 16293 Aa 


16:32:22.86 


-24:28:36.6 


980+200 


0.8x0.7 


0.05 


IRAS 16293 Ab 


16:32:22.88 


-24:28:36.3 


430+90 


0.9x0.6 


0.02 


IRAS 16293 B 


16:32:22.62 


-24:28:32.5 


3130+630 


0.6x0.4 


0.15 


L483 


18:17:29.94 


-04:39:39.3 


140+28 


2.2x0.8 


0.02 



^Center position, flux, and FWHM size of the dust continuum sources derived from Gaussian 
fitting. 

''Total gas mass derived from the SMA 850 /im dust continuum observations. Note that the gas 
mass is from compact emission only, not the whole envelope. 
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4. DISCUSSION 

4.1. Multiplicity in the Class Phase 

Before analyzing the binary/multiple frequency in the Class protostellar phase, it is 
important to note how binary/multiple systems are defined in this work. Strictly speak- 
ing, only systems that show bound motions can be referred to as binary/multiple systems. 
However, compared with well-studied MS/PMS binaries (for which the orbital motion can 
be determined by direct imaging or spectroscopic monitoring at optical and near-infrared 
wavelengths; see Duchene et al. 2007), these kinds of systems are extremely rare in the 
protostellar phase (see, e.g., Rodriguez 2004). The main reasons are: (1) the sample of 
protobinary systems is relatively small, (2) it is difficult to directly observe the protostellar 
phase, due to the limitation of angular resolution obtained at millimeter wavelengths, and 
(3) it is difficult to monitor bound motion in protobinary systems, as most known protobi- 
nary systems have large separations (100-1000 AU) and thus have long orbital periods. In 
this work, we use the following criterion to define binary/multiple systems: the separation 
or the closest separation (if there are more than two components in the systems) is less than 
5000 AU. We explain below why we choose 5000 AU as the limit. 

As introduced in Section 1, through initial fragmentation, large-scale (~ 0.1-1.0 pc) 
molecular clumps (or filaments) fragment into individual dense cores, which subsequently 
collapse to form individual stellar systems (i.e., fragmented components from initial clump 
fragmentation are not bound stellar systems but prestellar cores). The fragmented compo- 
nents from initial fragmentation generally have separations of more than 5000 AU, consistent 
with the typical Jeans length in large-scale molecular clumps/filaments (4000-5000 AU at T 
= lOK and nn^ = 10^-10^ cm~^). This is also in agreement with the results from the high 
angular resolution survey by Looney et al. (2000) that systems with separations < 6000 AU 
have common envelopes while wider systems have separated envelopes. Interestingly, by an- 
alyzing observational data in the Taurus region (e.g., Ghez et al. 1993; Gomez et al. 1993), 
Larson (1995) also found the existence of an 'intrinsic' scale (~ 0.04 pc or 8000 AU) in binary 
star formation process, and this scale is found to be consistent with the Jeans length in large 
molecular clouds. Therefore, we consider that components with separations > 5000 AU are 
likely not relevant for the formation of bound stellar systems. This is also the reason why 
we list several well-known sources, e.g., NGC 1333 IRAS2A & 2B (separation ~7500AU) 
and IRAS4A & 4B (separation ~7200AU), as separate targets in this work (see Table 1). 

We use the densities of YSOs and dense cores in the Perseus molecular cloud to esti- 
mate the likelihood of identifying a background or foreground source as a possible binary 
(or multiple) companion in this study. Many of the targets (14 of 33, or 42%) in this sam- 
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pie are located in the Perseus cloud, and most of the rest of the sources are in regions 
with lower protostellar densities. In general, main sequence and evolved stars foreground 
and background to the cloud have very small amounts of (or no) circumstellar dust, and 
should be undetectable by this interferometric submillimeter and millimeter continuum sur- 
vey. Therefore, the only possible "contaminating" sources in this study should be young 
stars with dusty circumstellar envelopes associated with the cloud. In Perseus, there are a 
total of 400 YSOs in an area of 3.86 deg^ (see J0rgensen et al. 2006). Assuming that these 
YSOs all have detectable dust continuum emission at the SMA, the density of potential dust 
continuum sources is 103.6 YSOs/deg^. We therefore estimate the probability of there being 
a contaminating YSO within a 5000 AU (21" at the distance of Perseus) radius that could 
mistakenly be identified as a companion to be 0.01. Another estimate of the contamination 
likehhood can be obtained using the density of dense cores detected in previous single-dish 
dust continuum surveys. Hatchell et al. (2005; SCUBA 850 /im) and Enoch et al. (2006; 
Bolocam 1.1mm) found 91 (within an area of 3deg^) and 122 (within an area of 7.5 deg^) 
dense cores in Perseus, respectively. From these two large-scale surveys, the probabilities are 
estimated to be 3.2 x 10~^ and 1.7 x 10^'^, respectively. Hence, it is very unlikely that any 
of the identified companions in this work are background or foreground sources. 

Nevertheless, we want to note that further high angular resolution kinematical data are 
important for complementing the separation criterion used here. For example, for relatively 
wide protobinary systems (> 1000 AU), molecular line observations are needed to verify if the 
components therein have similar systematic velocities, in order to insure that the components 
resolved therein are physically associatecTl For close systems (on the scale of 100 AU), 
high angular resolution line (e.g., CO or II2) observations are needed to distinguish if each 
component drives its own outflow, in order to be sure that each component is associated 
with a protostellar object (and not a dense, starless, condensation). 

In this work, we adopt the following terminology and definitions to describe the com- 
plexities of binary/multiple stars: 

multiplicity frequency (MF) = , (4) 

companion star fraction (CSF) = — — -, (5) 



■^For several wide protobinary systems (e.g., L1448N, L1448C, IRAM 04191, CG30, and L723), high 
angular resolution molecular line observations have already found that the components resolved therein have 
similar systematic velocities. 
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where S, B, T, and Q are the numbers of single, binary, triple, and quadruple systems, 
respectively. The 'multiplicity frequency', suggested by Reipurth & Zinnecker (1993), repre- 
sents the probability that any system is a binary/multiple system. The 'companion star frac- 
tion', on the other hand, looks at the individual stars and counts the number of companions 
relative to the sample size. It should be noted that many previous works used 'compan- 
ion star fraction' rather than 'multiplicity frequency' in their analysis (see, e.g., Haisch et 
al. 2004; Duchene et al. 2004; Connelley et al. 2008b). Actually, as noted by Reipurth & 
Zinnecker (1993), the most precise description of multiplicity consists of explicitly listing the 
number of singles, binaries, triples, and quadruples, like S:B:T:Q, although this description 
(as well as others) does not account for hierarchicality. 

In this SMA survey toward 33 protostars, 22 sources show signatures of binarity or 
multiplicity (see Appendix A for detailed descriptions of the individual sources). Table 6 
lists the names of these sources, as well as their basic observed parameters (separation and 
circumstellar mass ratio). If we adopt 5000 AU as the limit, source CB 17 (21'.'4 or 5350 AU) 
is beyond this limit and will not be regarded as a "bound" binary system in this work. 



Therefore, the numbers of singleq^ binarieqj triples ^'^, and quadruples'^ are considered to 
be 12, 14, 5, and 2, respectively, which is S:B:T:Q = 0.86:1.00:0.36:0.14 when normalized 
to the number of binaries. The overall multiplicity frequency (MF) and companion star 
fraction (CSF) derived in the sample, without correcting for completeness (see below), is 
0.64 ±0.08 and 0.91 ±0.05, respectively, with the separations in the range between 50 AU 
(which is the smallest separation we detect between two components) and 5000 AU (the limit 
set in this work). The statistical standard errors are estimated as [MF(1— MF)/N]-'^/^ and 
[CSF(1— CSF)/N]'/^ (where N is the total number of sources in the sample), respectively. 

It is important to note that the angular resolution is not uniform in this SMA survey 
(see Tables 2 & 3). The median angular resolution in the survey is 2'.'5, corresponding to 
a median linear resolution of 600 AU in nearby clouds (see Figure 1). There are very likely 
more binary /multiple systems with separations less than 600 AU unresolved in this survey. 
Therefore, this survey is incomplete across the observed separation range. However, it is 
difficult to correct for this incompleteness at present, because we would have to adopt as- 



8L1448 IRS2, NGC1333 IRAS2B, Blc, CB 17, L1521F-IRS, L1527IRS, HH114 MMS, B228, Oph 
MMS126, L483, B335, and L1157 

9L1448C, NGC1333 IRAS2A, IRAS4A, IRAS4B, IRAS 03282, Bib, HH 211 MMS, IC348MMS, 
IRAM 04191, 0MC3 MMS6, NGC2024-FIR5, HH25MMS, CG30, and L723VLA2 

10L1448N, PER 065, HH212MMS, VLA1623, and IRAS 16293 

"NGC 1333 SVS 13 and L1251B 



sumptions regarding to the separation distributions of Class protobinary systems, which 
are unfortunately unknown yet (see further discussion in §4.2). Consequently, the multi- 
plicity frequency and companion star fraction derived from this work could be considered as 
lower limits over the separation range observed here. 

Recently, Maury et al. (2010) argued that binary systems are relatively rare in the 
Class phase, as only one object shows a possible binary detection in their small PdBI 
survey toward five Class objects. However, two protobinary systems detected in the SMA 
survey presented here, L1448C (see Appendix A.l) and IRAM 04191 (see Appendix A. 8), 
are missing in their PdBI A-configuration (the most extended configuration) observations 
due to the low flux-ratio sensitivity (only 0.6-0.9; see Maury et al. 2010). Therefore, taking 
both L1448C and IRAM 04191 into account, the re-calculated multiplicity frequency in their 
sample is actually 0.60, which is consistent with the multiplicity frequency derived in this 
SMA survey (0.64). 



4.2. The Distribution of Separations 

Figure 19 shows the distribution of separations for the Class binary/multiple systems 
observed in this work (see Table 6). It must be noted that these separations, measured in 
the SMA (sub-) millimeter images, are the projected separations between the components in 
individual systems. However, as discussed by Reipurth & Zinnecker (1993), using reasonable 
statistical assumptions, one could convert the mean projected separation (/) to the mean 
semi-major axis (a) with the ratio of //a ~ 0.79-0.85 (see Reipurth & Zinnecker 1993 and 
references therein). We therefore consider that there is no significant error by assuming that 
the overall statistical distribution of observed projected binary separations is representative 
of the actual distribution of the semi-major axes. In this work, we do not use these ratios 
to convert observed projected separations to the 'semi-major' axes, and instead show the 
original results from the observations. 

As shown in the top of Figure 19, the distribution of separations for Class protostars 
shows a general trend in which the companion star fraction increases from wide separations to 
close separations. Quantitatively, the companion star fraction is 0.42 ±0.09 for separations 
from 50 AU to 1700 AU, 0.27 ± 0.08 from 1700 AU to 3350 AU, and 0.21 ± 0.07 from 3350 AU 
to 5000 AU. The median separation in this survey is 1800 AU. Previous large-scale surveys 
toward MS and PMS stars found that the distributions of separations could be well fitted 
with log-normal functions with peaks at less than 100 AU (see, e.g., Duquennoy & Mayor 
1991 and Patience et al. 2002). As shown in the bottom of Figure 19, the distribution 
of separations for Class protostars does not exhibit a log-normal like distribution, and 
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Table 6: Separations and mass ratios of binary/multiple systems observed in this work 



Source 


; Separation ["] 


Separation [AU] 


Mow/.Mhigh 


L1448N (A, B^ 


7.4±0.3 


1800 


0.21 ±0.06 


L1448N (A, NW^ 


16.5 ±0.3 


4000 


0.50 ±0.15 


L1448C (N, S^ 


7.0 ±0.5 


1700 


0.10±0.05 


PER 065 (A, B^ 


3.9±0.4 


940 


0.25 ±0.10 


PER 065 (B, C^ 


7.5 ±0.4 


1800 


0.73 ±0.22 


NGC1333IRAS2A (Al, K2] 


1.5±0.2 


360 


0.05 ±0.02 


NGC1333SVS13 (A, VLA3^ 


5.0±0.3 


1200 


0.05 ±0.02 


NGC1333SVS13 (VLA3, B^ 


10.4 ±0.3 


2500 


0.05 ±0.02 


NGC1333SVS13 (B, C^ 


19.8 ±0.3 


4750 


0.11 ±0.03 


NGC1333IRAS4A (Al, A2] 


1.4±0.2 


340 


0.50±0.15 


NGC1333IRAS4B/C 


10.5 ±0.2 


2500 


0.27±0.10 


IRAS03282+3035 (MMSl, MMS2^ 


1.6±0.2 


400 


0.10±0.05 


Per-Bl-b (N, S^ 


17.4 ±0.5 


4200 


0.71 ±0.21 


HH211 (SMMl, SMM2^ 


0.31 ±0.02 


74 


0.30±0.10 


IC348 (MMSl, MMS2^ 


9.8±0.4 


2350 


0.25 ±0.10 


'^CB 17 (IRS, MMS^ 


21.4±0.5 


5350 


0.57±0.17 


IRAM 04191+1522 (IRS, MMS^ 


8.0 ±0.4 


1100 


0.56±0.17 


OMC3MMS6 (SMMl, SMM2^ 


10.8 ±0.2 


4750 


0.41 ±0.12 


NGC2024 FIR 5 (SMMl, SMM2^ 


4.1±0.2 


1800 


0.39 ±0.12 


HH212 (SMMl, SMM2^ 


0.53 ±0.05 


230 


0.08 ±0.03 


HH212 (SMMl, SMM3^ 


1.2±0.1 


530 


0.10 ±0.03 


''HH25 (SMMl, SMM2^ 


12.7±0.3 


5600 


0.84 ±0.25 


HH25 (SMM2, SMMS^ 


10.9 ±0.3 


4800 


0.35 ±0.10 


CG30 (N, S^ 


21.8±0.5 


4400 


0.38±0.11 


VLA 1623 (W, E^ 


9.9±0.2 


1200 


0.20 ±0.06 


VLA 1623 (EA, EB^ 


1.1±0.1 


130 


1.0±0.3 


IRAS16293-2422 (A, B^ 


5.3±0.2 


640 


0.50 ±0.15 


IRAS16293-2422 (Aa, Ab^ 


0.42 ±0.05 


50 


0.44 ±0.13 


L723VLA2 (MMSl, MMS2^ 


3.5±0.3 


1050 


1.0±0.3 


L1251B (IRSl, SMA-N^ 


8.9 ±0.4 


2700 


0.35 ±0.11 


L1251B (SMA-N, SMA-S^ 


8.4±0.4 


2500 


0.61 ±0.18 


L1251B (SMA-S, IRS2^ 


15.6 ±0.4 


4700 


0.80 ±0.25 



° CB 17 is not regarded as a "bound" binary system in this work. 
^ The separation between sources SMM 1 and SMM 2 is beyond the Hmit of 5000 AU, and hence we regard 

HH25MMS as a binary system (SMM 2 and SMM 3). 
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rather shows a trend where most companions are found with separations around 1800 AU. 
We ascribe this distribution to the hmited angular resolution of the SMA observations (i.e., 
the completeness limit of the survey, see below), which allowed us to mostly detect binary 
separations in the order of 1000 AU. Further higher angular resolution large-scale surveys 
toward Class protostars are needed, in order to achieve a wider distribution of separations 
(from ~0.1AU to ~10^AU) for protobinary systems. 

To better understand the evolution of binary properties, it is important to compare the 
separation distribution of Class protobinary systems with those derived for Class I, PMS, 
and MS binaries. Yet, it must be noted in advance that MS and PMS binaries, as well as 
Class I binaries, have been well studied over the past two decades thanks to the availability 
of both large samples and high angular resolution obtained with large optical/near-infrared 
telescopes. The sample size and angular resolution in this SMA survey are an order of 
magnitude less than the optical and near-infrared surveys. Therefore, large uncertainties 
may remain in our comparisons due to the small sample size and limited angular resolution 
in this work. 

In a large-scale survey toward Class I YSOs, Connelley et al. (2008a, b) found the num- 
ber of single, binary, triple, and quadruple stars is S:B:T:Q = 122:51:12:4, with separations 
ranging from about 50 AU to 5000 AU. Based on the large-scale survey of Duquennoy & 
Mayor (1991), with a similar separations range, the corresponding number found in solar- 
type MS stars is S:B:T:Q = 131:32:1:0 (see Connelley et al. 2008b). These two studies, which 
summarize the binary properties for Class I YSOs and MS stars, will be mainly used here for 
comparison, as their separation ranges basically match the observed separation range in this 
SMA survey. From these two large-scale surveys, the multiplicity frequency and companion 
star fraction are derived to be 0.35 ±0.03 and 0.46 ±0.04 for Class I YSOs, and 0.20 ±0.03 
and 0.21 ±0.03 for MS stars, respectively. Despite the incompleteness in the SMA survey 
and the difference between the sample sizes, the multiplicity frequency (0.64 ±0.08) and 
companion star fraction (0.91 ±0.05) found in Class protostars are 1.8 (for MF) and 2.0 
(for CSF) times larger than the values found in Class! YSOs, and 3.2 (for MF) and 4.3 
(for CSF) times larger than the values found in MS stars (see Figure 20). Furthermore, the 
observed fraction of high order multiple systems to binary systems in Class protostars (7:14 
or 0.50 ±0.09) is also larger than the fractions found in Class I YSOs (0.31 ±0.07) and MS 



stars (0.03 ±0.01)^^ To estimate how different the distribution of separations for Class 
protostars is from those distributions for Class I YSOs and MS stars, we use the equation 
given by Brandeker et al. (2006) in their Appendix B2 to calculate the probability that the 



^^The estimated fraction from the whole sample of Duquennoy & Mayor (1991; from ^0.1 AU to ^ 10^ AU) 
is about 20% (see also Tokovinin 2004). 
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results in a given separation range from two studies are consistent with each other. Using this 
method, the probabihties that Class protostars have the same distribution of separations 
as Class I YSOs and MS stars are estimated to be 1.9 x 10"^ and 4.1 x 10""^^, respectively. 

The distributions of separations derived for PMS binaries are in general similar to those 
for Class I binaries (see, e.g., Connelley et al. 2008b). Therefore, we consider that the 
comparisons between Class and PMS binaries will likely hold similar results to those derived 
from the comparisons between Class and Class I binaries. Indeed, for example, in the survey 
toward PMS stars, Chez et al. (1997) found a companion star fraction of 0.21 ±0.04 with 
separations from 150 AU to 1800 AU; with the same separation range, the companion star 
fraction found by Reipurth & Zinnecker (1993) is 0.16 ±0.02. The companion star fraction 
derived from this survey with the same separation range is 0.39 ±0.08, approximately two 
times higher, comparable to the ratios estimated between the Class and Class I binaries. 

As we indicated in §4.1, this SMA survey is incomplete across the observed separa- 
tion range, and the derived multiplicity frequency and companion star fraction in this work 
should be considered as lower limits. Therefore, the ratios estimated above should also be 
considered as lower limits. Considering that the median separation in this SMA survey is 
1800 AU (larger than the lowest linear resolution in the survey of 1500 AU), the survey can be 
regarded to be complete for protobinary systems with separations from 1800 AU to 5000 AU. 
Since, to our knowledge, this work is the first large-scale observational survey of multiplicity 
in Class protostars, the distribution of protobinary systems with smaller separations is not 
yet known. Hence, it is difficult to correct for our incompleteness to derive an accurate multi- 
plicity frequency or companion star fraction for Class protobinary systems with separations 
between 50 AU and 1800 AU. Nevertheless, we may adopt a simple method to 'correct' for 
completeness in this work. For separations smaller than 1800 AU, we assume that the com- 
pleteness in a separation range is given by the averaged fraction of the sources in the total 
sample that have linear resolutions intermediate between the higher and lower limits of the 
range. For example, 55% of the sample has a resolution better than the median resolution of 
600 AU, then we assume that we only detect ^^^°+^'^o^° of the total number of multiple systems 
in the range between 600 AU and 1800 AU. Using this method, the multiplicity frequency 
and companion star fraction are derived to be 0.33 ±0.08 and 0.42 ±0.08 for systems with 
separations from 1800 AU to 5000 AU (100% completeness), 0.15 ±0.06 and 0.24 ±0.07 for 
systems with separations from 600 AU to 1800 AU (78% completeness), while 0.15 ± 0.06 and 
0.24 ±0.07 for systems with separations from 50 AU to 600 AU (32% completeness). After 
correcting for completeness, the overall multiplicity frequency and companion star fraction 
in this SMA survey are 1.0 ±0.1 and 1.5 ±0.2, respectively, over the separation range from 
50 AU to 5000 AU. Compared to the data derived from Class! YSOs survey (Connelley et al. 
2008a, b), these 'corrected' values, result in a difference, in both the multiplicity frequency 
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and the companion star fraction, of a factor of three between Class protostars and Class I 
YSOs, which is indeed larger than the factor of two estimated above. However, it must be 
noted that these 'corrected' values are uncertain. Further high angular resolution large-scale 
surveys toward Class protostars are needed, in order to draw a real statistically significant 
separation distribution for close protobinary systems. 

The decrease in both the multiplicity frequency and companion star fraction, as well 
as the fraction of high order multiple systems, suggests that binary properties change as 
stars evolve from the Class to Class I/PMS, all the way to the MS phases. As reviewed 
by Goodwin et al. (2007), the evolution of binary properties has been ascribed to two main 
mechanisms: (1) the rapid dynamical decay of young small- A^ systems (2 < A^ < 12; see, 
e.g., Reipurth 2000; Reipurth & Clarke 2001; Goodwin & Kroupa 2005); and (2) dynamical 
destruction of loose binary systems in dense clusters (see, e.g., Kroupa 1995). Here, we 
consider that small- A^ decay is likely the main process for the evolution of binary properties, 
because (1) most sources in this sample, as well as those in the Class I, PMS, and MS surveys 
discussed above, are not in dense cluster environments, and (2) the time scale for the decrease 
in binary fraction caused by destruction in dense cluster (a few 10 Myr) is much greater than 
the lifetimes of Class protostars (~ 10^-10^ yr) and Class I YSOs (~ 10^ yr). As estimated 
by Goodwin et al. (2007), a triple system is unstable to decay with a general half-life of a few 
10^ yr, which is consistent with the lifetime estimated for Class protostars. Thus, most of 
the decays/ejections are expected to happen during the Class phase (e.g., Reipurth 2000). 
This small-A^ decay will modify the binary properties by reducing the overall multiplicity 
frequency and companion star fraction, and reducing the fraction of high order multiple 
systems to binary systems. All have been evidenced by the comparisons between this work 
and previous large-scale surveys toward Class I and MS/PMS binaries. 

In short, we find that approximately two-thirds of Class protostars are binary or 
multiple systems in this SMA survey, and this could only be a lower limit to the fraction of 
Class protostars formed in such systems. The results derived from this work, in concert 
with the comparison between this survey and large surveys toward Class! and MS/PMS 
binaries, is consistent with the possibility that most, if not all, stars are formed in binary or 
multiple systems and that single stars result from the decay of multiple systems, as suggested 
by Larson (1972). 



4.3. The Distribution of Mass Ratios 

We have derived the circumstellar gas masses for the sources in the sample using the 
(sub-) millimeter dust continuum emission (see Sections and Tables 4 & 5). We note that 
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these values are uncertain and generally cannot be compared among the sources in the 
sample, because different sources were observed differently (i.e., with different angular res- 
olution, Mf -coverage, and different resulting sensitivity), which probe different size scales 
and amounts of material. Relative masses between components within individual systems 
are more reliable. Assuming, for simplicity, that sources in a multiple system have similar 
temperatures, then we can obtain the circumstellar mass ratio from the continuum emission 
flux ratio of the sources (see Equations 1-3). Table 6 lists the circumstellar mass ratios of the 
binary/multiple systems observed in this work, under the assumption that the circumstellar 
mass is proportional to the (sub) millimeter continuum flux. It must be noted that this 
circumstellar mass ratio does not necessarily reflect the mass ratio of the hydrostatic cores 
nor that of the final stars, but still provides a valuable comparison between the components 
in these systems. Figure 21 shows the distribution of circumstellar mass ratios for the pro- 
tobinary systems in this work. The distribution appears to be flat like that of more evolved 
long-period PMS and MS binary stars (see, e.g., Reipurth & Zinnecker 1993; Halbwachs et 
al. 2003). We find that 67% ±8% of protobinary systems have circumstellar mass ratios 
below 0.5, i.e., unequal masses systems are more common than equal masses systems in the 
sample. This result suggests that unequal-mass systems may be preferred in the process of 
(wide) binary star formation (separations on the order of 1000 AU). 

Figure 22 shows the distribution of circumstellar mass ratios versus separations. The 
distribution is somewhat scattered, similar to the distribution of PMS binaries (see, e.g., 
Reipurth & Zinnecker 1993). Numerical simulations predict that close binary systems (sepa- 
rations on the order of 1 AU) are likely to have mass ratios near unity (see Bate 2000), which 
is indirectly supported by statistical studies of MS binary systems (see, e.g., Halbwachs et 
al. 2003). However, we do not find this trend in Figure 22. It is probably because (1) the 
sample size in this work is still small, and (2) we cannot resolve close binary systems in 
which components may prefer similar masses. 



4.4. Sequential Fragmentation in Star Formation 

Different physical processes at different evolutionary stages may trigger different frag- 
mentation processes (e.g., Tohline 2002; Machida et al. 2008). At different evolutionary 
stages, the varying properties of a collapsing core (e.g., core temperature and density) will 
result in different fragmentation scales. As introduced in Section I, depending on the evo- 
lution of a collapsing core, fragmentation models can be divided into initial fragmentation, 
prompt (isothermal) fragmentation, adiabatic fragmentation and secondary fragmentation. 
Suggested by these models, the fragmentation of a core could take place sequentially during 
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the core's different collapsing phases, which would result in a hierarchical system eventually. 
However, it is yet unclear whether the fragmentation at one specific phase is dominant over 
other phases, and whether most collapsing cores go though all these fragmentation processes. 
Based on the extensive single-dish and interferometric observations, we show in Figure 23 a 
suggested sequential fragmentation picture. This is a rough outline and further efforts are 
needed to collect corresponding samples and to derive their kinematical properties, in order 
to understand the fragmentation mechanisms at different core's collapse phases in detail. 

The initial fragmentation of a large-scale clump or filament (~ 0.1-1.0 pc) occurs before 
the collapse of individual prestellar cores, and results in cores with separations of about 
10^-10*^ AU. This large-scale fragmentation has been frequently observed with single-dish 
millimeter telescopes (see, e.g., Kauffmann et al. 2008; Launhardt et al. 2010), and it 
is likely controlled by the combined effects of gravity and turbulence. Figure 23 shows a 
possible example of the initial clump fragmentation, CB 246 (see Launhardt et al. 2010 for 
more details), where the extended large-scale clump has fragmented into separated cores. 
These cores are invisible in the Spitzer Space Telescope (Spitzer) infrared images and may 
represent the prestellar cores which will further collapse and fragment. 

Most of the protobinary systems that have been observed with current millimeter ar- 
rays have separations in scales of 1000 AU (e.g., Looney et al. 2000, and this work). These 
systems appear to support the scenario of prompt fragmentation, which occurs at the end 
of the isothermal collapse phase (see Tohline 2002). Nevertheless, these observed systems 
are mostly made of Class objects that have already gone through the fragmentation pro- 
cess, and thus these systems do not exhibit the initial conditions of a fragmenting core. An 
observational example representing the prompt fragmentation may be the prestellar core 
RCrA SMM lA (see Figure 23), in which multiple faint condensations were discovered with 
separations between 1000 AU and 2000 AU. These condensations are invisible in the deep 
infrared images, have extremely low bolometric luminosities (< 0.1 L0) and temperatures 
(< 20 K), indicating that these are young sources have yet to form protostars, and therefore 
represent the earliest phase of core fragmentation observed (see Chen & Arce 2010 for more 
details). More recently, based on the kinematic data collected at different interferometric 
arrays, Chen et al. (2012a) found that most protostellar cores with binary systems (sep- 
arations scale ~ 1000 AU) formed therein have ratios of rotational energy to gravitational 
energy /^rot > 1%, which is consistent with theoretical simulations (see, e.g.. Boss 1999). 
This suggests that the level of rotational energy in a dense core plays an important role in 
the prompt fragmentation process. 

Unfortunately, there is a lack of observations of protobinary systems with separations 
between 10 AU and 100 AU, which could be used to study the adiabatic fragmentation see- 
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nario (see Machida et al. 2008). This can be easily explained by the fact that current 
millimeter interferometers (e.g., SMA and PdBI) normally reach 1-2" angular resolution 
under general conditions, and thus mostly resolve protobinary systems with separations of 
100-200 AU or larger in nearby star-forming clouds (e.g., Perseus and Ophuichus). In Fig- 
ure 23, we show the eSMA 850 /im image of source IRAS 16293A, which represents the best 
angular resolution that can be achieved at the SMA. In IRAS 16293 A, a close binary system 
with a separation of 50 AU is revealed (see also Appendix A. 19). We speculate that the close 
binary system in IRAS 16293A was formed through adiabatic fragmentation, which took 
place after the initial clump fragmentation (that formed the IRAS 16293 and IRAS 16293 E 
cores; see, e.g.. Figure 8 in J0rgensen et al. 2008) and the prompt (isothermal) fragmentation 
(that formed the IRAS 16293 A-B system; see Figure 15d). We note that the protobinary 
system VLA 1623 (see Appendix A. 17 and Figures 15a-b) may show another example of 
adiabatic fragmentation. Here, the VLA 1623 system appears to have fragmented from the 
Oph-A core through initial clump fragmentation and have later fragmented into VLA 1623 
West & East (through prompt fragmentation); VLA1623 East seems to have subsequently 
gone through adiabatic fragmentation to form VLA 1623 East A and B. As we discussed in 
Appendix A. 17, source VLA 1623 B could also be a binary system, which needs to be verified 
by further higher angular resolution observations. 

We believe that routine observations at angular resolution better than O'.'l will reveal 
more multiple protobinary systems with smaller separations, which will then allow us to study 
in detail a statistically significant number of protostellar binary /multipe systems with a wide 
range of separations (from 1 AU to 10^ AU). The Atacama Large Millimeter/submillimeter 
Array (ALMA), capable of attaining an angular resolution down to 10 mas (when fully com- 
pleted), will provide a breakthrough in our knowledge of binary star formation. It will 
certainly help us find closer protobinary systems with separations in the scale of 1-10 AU, 
and will even allow us to test the secondary fragmentation scenario (see, e.g., Machida et al. 
2008), by which stellar systems with separations of 0.01-0.1 AU can be formed. 



5. SUMMARY 

We present SMA 1.3 mm and 850 /im dust continuum data toward a sample of 33 Class 
protostars in nearby clouds (distance < 500 pc), which thus far represents the largest survey 
toward protostellar binary/multiple systems. The median angular resolution in the survey 
is 2'.'5, while the median linear resolution is approximately 600 AU. The main results of this 
work are summarized below. 

(1) Compact millimeter and/or submillimeter dust continuum emission is detected from 
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all sources in the sample, which likely originates from the inner envelope and/or circumstellar 
disk. The measured 1.3 mm fluxes in the sample range from 1.1 mJy to 3.0 Jy, with a median 
value of 100 mJy, while the 850 /xm fluxes range from 20mJy to 3.2 Jy, with a median value 
of 430 mJy. Assuming that the 1.3 mm and 850 /xm dust continuum emission is optically thin 
and a dust temperature of 30 K for all sources, the estimates of the total gas masses are in 
the ranges of 0.001-2.4 Mq at A 1.3 mm and of 0.005-1.7 Mq at A850/im. 

(2) Twenty-one sources in the sample show signatures of binarity/multiplicity, with 
separations ranging from 50 AU to 5000 AU. Of them, four sources are newly-discovered 
binary/muhiple candidates (PER 065, NCG 1333 IRAS2A, IC348 MMS, and HH25SMM). 
Quantitatively, the numbers of singles, binaries, triples, and quadruples in the sample are 
12, 14, 5, and 2, respectively. The estimated multiplicity frequency and companion star 
fraction are 0.64 ± 0.08 and 0.91 ± 0.05, respectively, with no correction for completeness. We 
consider these values as lower limits, as it is likely that there are unresolved binary/multiple 
systems with small (< 600 AU) separations in this survey. 

(3) The derived multiplicity frequency and companion star fraction for Class proto- 
stars are approximately two times higher than those values found in the binary surveys 
toward Class! YSOs and pre-main sequence stars, and approximately three (for multiplicity 
frequency) and four (for companion star fraction) times larger than the values found among 
main sequence stars, with a similar range of separations. Furthermore, the observed fraction 
of high order multiple systems to binary systems in Class protostars (0.50 ±0.09) is also 
larger than the fractions found in Class I YSOs (0.31 ± 0.07) and main sequence stars (< 0.2). 
The decrease in both the multiplicity frequency/companion star fraction and the fraction 
of high order multiple systems clearly suggests that binary properties evolve as protostars 
evolve, as predicted by numerical simulations. 

(4) The distribution of separations for Class protostars shows a general trend in which 
the companion star fraction increases with decreasing companion separation. Quantitatively, 
the companion star fraction is 0.42 ± 0.09 for separations from 50 AU to 1700 AU, 0.27 ± 0.08 
from 1700 AU to 3350 AU, and 0.21 ±0.07 from 3350 AU to 5000 AU, with no correction for 
completeness. The median separation in this survey is 1800 AU (larger than the lowest 
linear resolution), and thus the survey is regarded to be complete for protobinary systems 
with separations from 1800 AU from 5000 AU (over this separation range the multiplicity 
frequency and companion star fraction are 0.33 ±0.08 and 0.42 ±0.08, respectively). After 
correcting for incompleteness for separations less than 1800 AU with a simple method, the 
overall multiplicity frequency and companion star fraction in this survey are 1.0 ±0.1 and 
1.5 ±0.2, respectively, over the separation range from 50 AU to 5000 AU. 

(5) The distribution of circumstellar mass ratios for the protobinary systems in this 
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survey appears to be flat, like that of more evolved long-period main-sequence and pre-main 
sequence binary stars. We find that 67% ± 8% of the protobinary systems have circumstellar 
mass ratios below 0.5 (i.e., unequal masses are much more common than equal masses). 
This implies that unequal-mass systems are preferred in the process of binary star formation 
(with separations in the scale of 1000 AU). 

(6) The high angular resolution survey in this work, in concert with previous extensive 
single-dish observations, strongly support that the fragmentation of molecular cloud cores, 
at different core collapse phases, is the main mechanism for the formation of binary/multiple 
stars. Based on these observations, we suggest an empirical sequential fragmentation picture 
for binary star formation. 

We thank the anonymous referee for providing many insightful suggestions and com- 
ments, which helped us improve this work greatly. We thank the SMA staff for technical 
support during the observations and their maintenance of the SMA archival data. X.C ac- 
knowledges the support of the Thousand Young Talents Program of China. This material is 
based on work supported by NSF grant AST-0845619 to H.G.A. JEP has received funding 
from the European Community's Seventh Framework Programme (/FP7/2007-2013/) under 
grant agreement No 229517. 
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A. DESCRIPTION OF INDIVIDUAL SOURCES 

In this Appendix, we provide a brief description of the individual sources, and discuss 
their SMA results in more detail. We note that, for sources L1448C, NGC 1333 IRAS2A & 
2B, NGC 1333 IRAS4A & 4B, L1527, L483, B335, and L1157, the SMA results (dust con- 
tinuum and molecular lines) were presented by J0rgensen et al. (2007). The dust continuum 
results (images and fluxes) derived independently in this work are comparable with those in 
J0rgensen et al. (2007). We refer readers to that paper for more details on the line data 
of these sources. Other sources for which their SMA data have been published elsewhere 
include: HH211 MMS (Lee et al. 2007b; 2009; 2010), NGC2024-FIR5 (Alves et al. 2011), 
HH212 MMS (Lee et al. 2006; 2007a; 2008), L723 VLA2 (Girart et al. 2009), and L1251B 
(J. Lee et al. 2006; 2007). 



A.l. L1448 Region 

L1448 is a dark globule located on the west end of the Perseus cloud, which has been 
observed extensively in the past two decades (see, e.g.. Bally et al. 2008, and references 
therein). Figure 3a shows the SCUBA 850 ;um image of the L1448 region. The region 
contains three well-known Class protostars: L1448 IRS2, L1448N, and L1448C (see, e.g., 
Barsony et al. 1998; O'Linger et al. 2006). 

L1448 IRS2 was suggested to be a protobinary candidate (with a 10" separation) by 
Wolf-Chase et al. (2000) based on the NRA0-12m CO (1-0) observations that mapped two 
distinct molecular outflows from IRS2. Volgenau et al. (2002) also claim the detection of a 
binary system in IRS2 using BIMA, but no results (e.g., spatial separation or image) have 
been presented yet. However, the SMA 1.3 mm dust continuum observations with an angular 
resolution of 3'.'4 x 2'.'5 (see Figure 3c), as well as Spitzer infrared observations (e.g., Chen 
et al. 2010), do not flnd any companions around IRS2 (see also O'Linger et al. 2006). We 
therefore conclude that L1448 IRS2 is a single protostar. 

L1448N is also referred to as L1448IRS3 by several authors. However, we prefer 
using the name L1448N, because the IRAS source IRS3 in fact consists of two protostars, 
L1448N and L1448C (see O'Linger et al. 2006). L1448N is the brightest source at infrared 
and millimeter wavelengths in L1448, and has been resolved into three distinct sources by 
previous interferometric observations (e.g., Looney et al. 2000). The three sources are 
named L1448NA, L1448NB, and L1448NW (O'Linger et al. 2006). The SMA 1.3 mm 
dust continuum image of L1448N shows these three distinct sources (see Figure 3b). The 
dust continuum emission is dominated by source NB, which is located 7'.'4 ± 0'.'3 south of 
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source NA. Source L1448NW, located approximately 20" northwest of the NA-NB pair, has 
relatively weaker dust continuum emission in the interferometric observations. 

L1448C, also known as L1448-mm, is one of the best-studied Class protostars. L1448C 
is well-known for driving a young, highly collimated, and high- velocity bipolar outflow (e.g., 
Bachiller et al. 1990; Hirano et al. 2010). Spitzer infrared observations resolved L1448C 
into two separated components, with an angular separation of about 8" (see J0rgensen et 
al. 2006). The SMA observations at both 850 /xm and 1.3 mm continuum bands detect these 
two components (see Figures 3d & 3e). The dust continuum emission is dominated by the 
northern component, and only weak dust continuum emission is found from the southern 
component (5 a detection). The northern component is the driving source of the well-known 
high- velocity bipolar outflow (J0rgensen et al. 2007; Hirano et al. 2010). Maury et al. 
(2010) argued that the southern component is not a protostar but a dusty clump associated 
with the outflow driven by the northern component. However, we note that a relatively low- 
velocity bipolar outflow associated with the southern component has been found in the SMA 
CO (2-1) observations (see J0rgensen et al. 2007) and CO (3-2) observations (see Hirano et 
al. 2010), which clearly indicates that the southern component is a protostellar object. 



A.2. PER065 

The protostellar core PER 065 was discovered by the dust continuum survey in the 
Perseus cloud, and was classified as a Class protostar by Hatchell et al. (2007) and Enoch 
et al. (2009). Because of its low luminosity (0.03-0.22^0), it was also cataloged as one of 
the very low luminosity objects (VeLLOs) in Perseus (No. 065; see Dunham et al. 2008). 
The SMA 1.3 mm dust continuum observations resolve, for the first time, the faint core into 
three components, named A, B, and C, respectively (see Figure 4). The separations between 
components A and B and between B and C are 3'.'9 ± 0'.'4 and 7'.'5 ± 0'.'4, respectively. 



A.3. NGC 1333 Region 

NCG 1333 is a well-studied protostellar cluster located in the western part of the Perseus 
cloud (see, e.g., a review by Walawender et al. 2008 and references therein, and Arce et al. 
2010). Submillimeter and millimeter dust continuum observations have revealed several 
protostellar cores within the region (see, e.g., Sandell & Knee 2001), of which three cores, 
SVS 13, IRAS 2, and IRAS 4 are the most well-known and studied (see Figure 5a). 

Figure 5b shows the SMA 1.3 mm continuum image of SVS 13, in which three distinct 
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sources, referred to as A, B, and C, are found. The SMA results are consistent with the 
results found in other high angular resolution observations (e.g., Bachiller et al. 1998; Looney 
et al. 2000; Chen et al. 2009). Another weak continuum source, which is spatially coincident 
with the radio source VLA3 (Rodriguez et al. 1997) and the 3 mm dust continuum source 
detected by Looney et al. (2000) and Chen et al. (2009), is found in the SMA 1.3mm dust 
continuum image (see Figure 5c). Also note that, with the VLA centimeter observations, 
Anglada et al. (2000; 2004) resolved source A into a close binary system with an angular 
separation of 0'.'3. 

The SMA continuum images of IRAS 2 and IRAS 4 are shown in Figures 5d & 5e and 
Figures 5f & 5g, respectively. The continuum results of the two sources derived independently 
in this work are comparable with those in J0rgensen et al. (2007). We therefore refer readers 
to that work for more details about the SMA results. We note that IRAS 2A is resolved into 
two separated components in the high angular resolution 850 /im image (see Figure 5e), with 
a faint companion 1'.'5 ± 0'.'2 to the north of the main dust peak. In the CO observations, 
a quadrupolar outflow centered at IRAS2A was detected (see, e.g., J0rgensen et al. 2007). 
This quadrupolar outflow is likely two different outflows driven by the two binary components 
in IRAS 2A. Further high angular resolution observations are needed to confirm the tentative 
detection of the faint secondary source in IRAS 2A and to distinguish which source drives 
which outflow. 



A.4. IRAS 03282+3035 

IRAS 03282+3035 (IRAS 03282) is a young Class protostar located in the western part 
of the Perseus cloud, which has been studied by various groups using different molecular line 
transitions (see, e.g., Bachiller et al. 1991, 1994; Arce & Sargent 2006; Chen et al. 2007; 
Tobin et al. 2011). 

The OVRO 1.3 mm dust continuum observations (at 0'.'9 x 0'.'7 resolution) revealed two 
millimeter sources in IRAS 03282, with an angular separation of 1'.'5 (see Launhardt 2004). 
Although the SMA 1.3 mm observations (at 2'.'5 x 2'.'2) do not resolve the two components, 
the SMA 850 ;um observations (at 0'.'9 x 0'.'8) show two separated components (see Figure 6). 
We refer to the main continuum source as MMS 1 and the faint source as MMS2. In the 
SMA images, the angular separation between the two sources is 1'.'6 ± 0'.'2, similar to the 
result found in the OVRO 1.3 mm images. 
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A. 5. Perseus Bl-c and Bl-b 

Barnard 1 (Bl) is a dark cloud in the western part of Perseus. Single-dish dust contin- 
uum observations have revealed four protostellar cores in this cloud (see, e.g., Matthews & 
Wilson 2002, Hatchell et al. 2005, and Enoch et al. 2006), which are commonly referred to 
as Bl-a, -b, -c, and -d (see Figure 7a). 

The Class protostar Bl-c was studied with BIMA by Matthews et al. (2006). The 
BIMA 3.3 mm dust continuum observations (at 6'.'4 x 4'.'9 resolution) revealed a protostellar 
envelope with radius of ~2400AU and mass of 2.1-2.9 Mq. The SMA 1.3mm dust contin- 
uum observations, with higher angular resolution (3'.'4 x 2'.'9), resolve an inner envelope with 
radius of 700 ± 100 AU and mass of 0.24 ± 0.05 Mq (see Figure 7b). 

The Bl-b protostellar core is elongated in the north-south direction in the SCUBA 
images (see Figure 7a). The 3 mm dust continuum observations at the Nobeyama Millimeter 
Array (NMA) found two separate sources in the Bl-b core, referred to as Bl-bN and Bl-bS 
(Hirano et al. 1999). The SMA 1.3mm observations show similar results to those found 
by Hirano et al. (1999): the two sources are separated by 17'.'4±0'.'5 and have similar 
circumstellar masses (mass ratio 0.7 ±0.2). 



A.6. IC 348 MMS and HH 211 MMS 

IC 348 MMS is a Class protostar discovered in the southwestern part of the IC 348 
cluster (see Eisloffel et al. 2003). It drives a collimated molecular outflow, traced by II2 
emission (Eisloffel et al. 2003) and CO emission (Tafalla et al. 2006). The SMA 1.3mm 
dust continuum images of IC 348 MMS reveal, for the first time, two sources embedded in the 
common envelope detected by single-dish dust continuum observations (see Figure 8b). The 
angular separation between the two sources is 9'.'8±0'.'4, and the mass ratio is 0.25 ±0.10. 
The main continuum source MMS 1 is the driving source of the collimated molecular outflow. 

The Class protostar HH211 MMS is well-known for driving a highly collimated 
bipolar outflow, which has a narrow bipolar jet-like component in the high- velocity CO/SiO 
emission and a relatively wide-angle bipolar cavity in the low- velocity CO emission (see, e.g., 
Gueth & Guilloteau 1999; Hirano et al. 2006; Lee et al. 2007b). 

Based on very high angular resolution SMA observations (at 0'.'20 x 0'.'15), Lee et al. 
(2009) found that HH211 MMS actually contains two dust continuum sources, SMM 1 and 
SMM 2, with an angular separation of 0'.'31 ± 0'.'02 (see Figure 8c). Source SMM 1 appears to 
be the driving source of the collimated jet/outflow, while the nature of source SMM 2, which 
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has a mass of only 1.5-4 Mjup, is still uncertain. More interestingly, based on the reflection- 
symmetric wiggle of the HH 211 jet, Lee et al. (2010) suggested that source SMM 1 itself could 
be a very low mass protobinary with a separation of ~ 4.6 AU (the putative 4.6 AU companion 
was not directly detected, but inferred from the morphology of the outflow). Further higher 
angular resolution observations are needed to conflrm this interesting tentative detection. 
More details about HH211 MMS and its surrounding envelope can be found in Gueth & 
Guilloteau (1999), Lee et al. (2007b; 2009; 2010), and Tanner & Arce (2011). 



A.7. CB17 

CB 17 (also known as L1389) is a small, isolated, and slightly cometary-shaped dark 
cloud, which was classifled as a Bok globule by Clemens & Barvainis (1988). CB 17 has been 
studied by various groups using different molecular line transitions (e.g.,Pavlyuchenkov et 
al. 2006) and multi-wavelenghth dust continuum emission (e.g., Launhardt et al. 2010). 

The SMA 1.3 mm dust continuum observations reveal two sources within CB 17, which 
are separated by 21'.'4±0'.'5 in the northwest-southeast direction (see Figure 9a). The north- 
western dust continuum source is associated with an infrared source seen in the Spitzer 
images and is referred to as CB 17 IRS (see Chen et al. 2012b). No compact infrared emis- 
sion is detected from the southeastern source in the Spitzer bands from 3.6 to 70 /im, and 
we refer to this source as CB 17 MMS. Interestingly, the SMA CO (2-1) observations suggest 
that source CB 17 MMS, an object with an extremely low bolometric luminosity (< 0.04 Lq), 
is driving a low velocity (~2.5kms^^) molecular outflow. These characteristics have lent 
this source to be classifled as a candidate of the so-called flrst hydrostatic cores (see Chen 
et al. 2012b for more details). 



A.8. IRAM 04191+1522 

IRAM 04191+1522 (IRAM 04191) is a well-studied Class protostar located in the 
southern part of the Taurus molecular cloud (Andre et al. 1999). IRAM 04191 drives a 
large bipolar molecular outflow, features bright molecular line and dust continuum emission, 
and shows evidences of supersonic infall, fast rotation, depletion, and deuteration (see, e.g., 
Andre et al. 1999; Belloche et al. 2002; Belloche & Andre 2004; Lee et al. 2005). 

The SMA 1.3 mm dust continuum observations reveal two distinct continuum sources 
with an angular separation of 8'.'0 ± 0'.'4, embedded within the elongated core of IRAM 04191 
(see Figure 9b). The discovery of this binary system was reported by Chen et al. (2012a). 
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The southeastern source, associated with an infrared source seen in the Spitzer images 
(see Dunham et al. 2006), represents the well-known Class protostar and is referred to as 
IRAM 04191 IRS here. In contrast, the northwestern source has no infrared emission detected 
in the Spitzer images, and is referred to as IRAM 04191 MMS (see Chen et al. 2012a for 
more details). 



A.9. L1521F-IRS and L1527-IRS 

L1521F is a well-studied dense core in Taurus. This core shows high central density, infall 
asymmetry, molecular depletion, and enhanced deuterium fractionation, and was suggested 
as a highly evolved starless core by Crapsi et al. (2004). However, Spitzer observations 
found that L1521F harbors a very low luminosity protostar, refereed to as L1521F-IRS 
(see Bourke et al. 2006), which is in the early stage of the Class phase (Shinnaga et al. 
2009; Terebey et al. 2009). The SMA 1.3mm dust continuum observations of L1521F-IRS 
show a faint source (see Figure 9c). More than 95% of the flux around L1521F-IRS is resolved 
out by the SMA observations, compared to the flux detected in the IRAM-30m 1.2 mm maps 
(600±150mJy; Bourke et al. 2006). We also note that L1521F-IRS was observed at higher 
angular resolution with the IRAM-PdBI by Maury et al. (2010). In the PdBI 1.3 mm dust 
continuum images, there is also only one object detected. 

L1527-IRS is another protostar in Taurus, which is associated with a far-infrared source 
and is probably more evolved than L1521F-IRS (see, e.g., Terebey et al. 2009). Based on the 
JCMT 800 /im continuum observations. Fuller et al. (1996) suggested L1527-IRS is a binary 
system with the secondary at about 20" from the main protostar. However, no evidence is 
found of this secondary source in the SMA 1.3mm dust continuum image (see Figure 9d). 
We therefore conclude that L1527-IRS is a single protostar (see also Maury et al. 2010). 



A.IO. HH 114 MMS 

HH 114 MMS is a Class protostar located in the western part of the A Orionis molecular 
shell, which drives a highly collimated bipolar outflow (see Arce & Sargent 2006). In the SMA 
1.3mm and 880 /xm dust continuum images (see Figure 10), source HH114 MMS remains 
as a single object at the angular resolution of 2'.'2 x 1'.'9, which is consistent with the OVRO 
2.7mm and 3.4 mm continuum observations (Arce & Sargent 2006). 



34 



A.ll. OMC3MMS6 

The MMS 6 core is the brightest source in submilhmeter and milhmeter wavelengths 
in the Orion molecular cloud-3 (OMC-3) region (Chini et al. 1997; Lis et al. 1998; John- 
stone & Bally 1999). The SMA 850 /im continuum images show two sources in the MMS 6 
core, referred to as SMM 1 and SMM 2, with an angular separation of l(y.'8 ± Qf!2 (see Fig- 
ure 11a). The northern continuum source SMM 1 is spatially coincident with a continuum 
source detected in early interferometric observations (see Matthews et al. 2005; Takahashi 
et al. 2009). The southern continuum source SMM 2 is associated with a continuum source 
detected in the BIMA 1.3mm observations by Matthews et al. (2005), but it was not de- 
tected in the multi- wavelength observations by Takahashi et al. (2009). We note that source 
MMS 6 was also observed in the SMA 0.9 mm continuum by Takahashi et al. (2009), but 
the sensitivity in their image was about 2 times worse than that in the SMA 850 /xm image 
presented here. 



A.12. NGC2024FIR5 

FIR 5 is a Class protostar located in the NGC 2024 cloud, the most active star- forming 
region in the Orion B cloud (see, e.g., a review by Meyer et al. 2008 and references therein). 
With a high velocity, collimated molecular outflow extending over 5' south of the core (Richer 
et al. 1992), FIR 5 is the brightest and also probably most evolved protostellar core among 
the dense cores in NGC 2024 (Mezger et al. 1992). 

In the SMA 850 /xm continuum images, the FIR 5 core is resolved into two main com- 
ponents, referred to as SMM 1 and SMM 2, with an angular separation of 4'.'1 ± 0'.'2 (see 
Figure lib). The SMA results are consistent with the results derived from the PdBI 3mm 
continuum observations, in which the FIR 5 core was also resolved into two main components 
(see Wiesemeyer et al. 1997). Nevertheless, in the BIMA 1.3mm deep observations, Lai et 
al. (2002) also found several weaker dust condensations around the two main components, 
which are not detected in the SMA images. We note that the SMA results of FIR 5 (includ- 
ing also dust polarization and CO emission results) were independently published by Alves 
et al. (2011), and similar dust continuum results were derived there. 



A.13. HH 212 MMS 

HH212 is a remarkable jet in the Orion L1630 cloud, discovered by Zinnecker et al. 
^1998). The jet is associated with a bipolar molecular outflow seen in CO lines (see, e.g.. 
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Lee et al. 2000; 2006; 2007a) and SiO lines (see, e.g., Codella et al. 2007; Lee et al. 2008). 
The driving source of the jet/outflow, referred to as HH 212 MMS, is a low-luminosity Class 
protostar associated with a cold IRAS source (IRAS 05413-0104; Zinnecker et al 1992). 

The SMA 1.3 mm and 850 /im dust continuum images of HH212 MMS were shown by 
Lee et al. (2006; 2007a; 2008). By re-analyzing the high angular resolution 850 /im images 
in Lee et al. (2008; angular resolution 0'.'35 x 0'.'32), a faint continuum component is found 
to the west of source HH 212 MMS (see Figure 12c). Hereafter, we refer to the main dust 
continuum component as MMS 1 and to the faint continuum component as MMS 2. The 
separation between the two components is measured to be 0'.'53±0'.'05. We also retrieved 
the SMA 230 GHz data of HH 212 MMS obtained in the Extended and Very Extended config- 
urations (the 230 GHz data published by Lee et al. 2006 were obtained only in the Compact 
configuration). The high angular resolution 1.3 mm dust continuum images (angular resolu- 
tion 0'.'8 X 0'.'7) show that source HH 212 MMS is extended in the east-west direction and a 
component is tentatively detected in the 1.3 mm images, which is spatially associated with 
the new component MMS 2 seen in the 850 /im images (see Figure 12d). 

We also note that another dust continuum source was tentatively detected in the PdBI 
1.4 mm continuum observations, which is about 1'.'7 to the southeast of the main dust con- 
tinuum source MMS 1 (see Codella et al. 2007). In the SMA 850 /im images, a weak dust 
continuum source (~7(T detection) is seen near this PdBI 1.4 mm source but 0'.'6 closer to 
the main continuum source MMS 1. We referred to this weak continuum source as MMS 3 
(see Figure 12c). Altogether, there may be three dust continuum sources found within the 
elongated envelope of HH 212 MMS. Further high angular observations are needed to confirm 
this tentative multiple system detection. 



A.14. HH25MMS 

The protostellar jet HH 25 is located in the northern part of the Orion B dark cloud. 
The jet is associated with a compact but coUimated molecular outflow (e.g., Gibb & Davis 
1998; Gibb et al. 2004). The driving source of the jet/outflow, referred to as HH25MMS, 
was revealed by a VLA 3.6 cm survey in this region (Bontemps et al. 1995), and was classified 
as a Class protostar with a luminosity of ~6Lq (Gibb & Davis 1998). 

In the single-dish dust continuum maps (IRAM-30m, Lis et al. 1999; JCMT/SCUBA, 
Phillips et al. 2001), HH25MMS is elongated in the north-south direction, with the VLA 
source located at the northern tip of the extension. In the SMA 870 /im dust continuum 
image, three distinct, previously unknown sources are revealed within the large-scale elon- 
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gated core, approximately aligned in the north-south direction (see Figure 13a). The three 
continuum sources, from north to south, are referred to as SMM 1, SMM 2, and SMM 3. The 
angular separations between sources SMM 1 and SMM 2 and between sources SMM 2 and 
SMM 3 are measured to be 12'/7 ± 0'.'3 and ICy/Q ± (y.'3, respectively. Source SMM 1 is spatially 
coincident with the VLA 3.6 cm source, and is the driving source of the HH25 jet/outflow. 



A.15. CG30 

CG 30 (also known as BHR 12; Bourke et al. 1995) is a bright-rimmed cometary globule 
located in the Gum Nebula region. In the SCUBA submillimeter (Henning et al. 2001) and 
ATCA 3 mm (Chen et al. 2008a) continuum observations, this globule was resolved into a 
wide protobinary system with an angular separation of ~ 22". 

The SMA 1.3 mm dust continuum image of CG30 shows two compact sources (see 
Figure 13b). The angular separation between the two sources is measured to be 21'.'8 ± O'.'S, 
consistent with the ATCA result. Following Chen et al. (2008a), we refer to the northern 
source as CG 30N (which is spatially coincident with IRAS 08076—3556) and to the southern 
source as CG30S. The two sources each drive their own bipolar protostellar jets/outflows, 
which are roughly perpendicular to each other (see Chen et al. 2008a, b). Further discussion 
of the gas kinematics (from ATCA N2H+ and SMA CO observations) and infrared emission 
(from Spitzer observations) in CG30 can be found in Chen et al. (2008a, b). 



A.16. B228 

B228 is a Class protostar located in the Lupus I cloud (Shirley et al. 2000). It drives 
a compact bipolar molecular outflow that extends in the northeast-southwest direction (see 
van Kempen et al. 2009). In the SMA dust continuum images, an elongated source is seen at 
both 1.3 mm and 1.1 mm (see Figure 14). The source is extended in the southwest direction, 
roughly following the direction of the blueshifted outflow lobe, suggesting that the outflow 
has a strong impact on the dust envelope, similar to L1157 (see, e.g., Gueth et al. 2003). 



A.17. VLA 1623 

VLA 1623 is the prototypical Class protostar (Andre et al. 1993), which drives a 
large scale bipolar outflow (Dent et al. 1995; Yu & Chernin 1997). With VLA centimeter 
continuum observations, Bontemps & Andre (1997) found a series of emission clumps, named 
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A, B, and C At high angular resolution, clump A was further resolved into two components 
(Bontemps & Andre 1997; Looney et al. 2000). 

In the SMA high angular resolution images (resolution &!6 x 0'/3; observed in 2007), 
three distinct dust continuum sources are detected (see Figure 15a). The three sources, 
aligning roughly from east to west, are referred to as VLA 1623 A, VLA 1623 B, and VLA 1623 
West, respectively. Sources A and B, the east pair with an angular separation of I'.'l iO'.'l 
(see Figure 15b), are spatially coincident with the BIMA 2.7 mm sources A and B in Looney 
et al. (2000) and the two VLA centimeter sources in clump A in Bontemps & Andre (1997). 
The mass ratio between sources A and B is close to 1.0 ± 0.3, similar to the result found by 
Looney et al. (2000). Source VLA 1623 West, located approximately 10" to the east of the 
A-B pair, is spatially coincident with clump B in Bontemps & Andre (1997). 

More recently, Maury et al. (2012) presented high angular resolution SMA data of 
VLA 1623 (observed in 2009) and argued that sources VLA 1623 B and VLA 1623 West are 
two dusty knots due to shocks along the outflow cavity wall of source VLA 1623 A (although 
they did not show any evidence that VLA 1623 A is the driving source of the VLA 1623 
outflow). However, in the Spitzer mid-infrared images, source VLA 1623 West is associated 
with a point-like infrared source, strongly suggesting that it is a young stellar object (X. Chen 
et al. in preparation). Furthermore, the SMA CO (2-1) data show a quadrupolar outflow 
centered at source VLA 1623 B, indicating that this source is not only a protostellar object 
but also a potential close binary system unresolved in the current SMA dust continuum 
observations (X. Chen et al., in preparation). 



A.18. Oph-MMS 126 (IRAS 16253-2429) 

Oph-MMS 126 is a Class protostar discovered in the SEST 1.2 mm continuum survey 
toward the p Ophiuchi molecular cloud, which is associated with a faint, cold IRAS source 
(IRAS 16253-2429) and drives a low-velocity bipolar molecular outflow (see Stanke et al. 
2006). In the SMA images, a faint dust continuum source is seen (about 5 a detection; see 
Figure 15c), which is associated with an infrared source seen in the Spitzer images (see 
Barsony et al. 2010). The total flux was measured to be ll±2mJy in the SMA 1.3mm 
dust continuum images. Compared with the total flux obtained in the SEST observations 
(411 mJy; Stanke et al. 2006), more than 97% of the dust continuum emission flux from the 
large envelope was resolved out by the SMA. 
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A.19. IRAS 16293-2422 

IRAS 16293-2422 (IRAS 16293) is a well-studied Class protostar located in the p Ophi- 
uchi cloud. Millimeter interferometric observations resolved this source into a binary system, 
referred to as IRAS 16293 A and B, with a separation of 5'.'2 (see, e.g., Mundy et al. 1992; 
Loony et al. 2000). In the high angular resolution VLA radio continuum observations, source 
A is further resolved into a chain of radio components (see Loinard et al. 2007a; Pech et al. 
2010), which are suggested to be related to a jet from source A. Recent ALMA Science ver- 
ification data revealed an inverse P-Cygni profile towards source B, which is direct evidence 
of infalling material close to the central protostar (Pineda et al. 2012). 

Figure 15d shows the SMA 1.3 mm continuum image of IRAS 16293, in which two con- 
tinuum sources, A and B, are detected with an angular separation of 5'.'3 ± 0'.'2. In the higher 
angular resolution eSMA 850 /xm images, source A is further resolved into two close con- 
tinuum sources (see Figure 15e). The separation between the two sources is 0'.'42 ± O'.'OS. 
Following Chandler et al. (2005), the two sources are referred to as Aa and Ab. 



A.20. L483 

L483 is an isolated dense core, associated with a strong infrared source (IRAS 18148—0440). 
It was classified as a Class O/I transition object by Tafalla et al. (2000). High angular reso- 
lution molecular line observations found a cold, quiescent core in L483 (see J0rgensen 2004). 
In the SMA 850 /im dust continuum observations, source L483 remains as single object at 
the angular resolution of 2'.'6 x 1'.'9 (see Figure 16), which is consistent with previous inter- 
ferometric millimeter continuum observations (see, e.g., J0rgensen 2004). 



A.21. L723VLA2 

L723 is a small isolated dark cloud. At the center of the cloud core is an IRAS point 
source (IRAS 19156+1906), which was classified as a Class protostar with a bolometric 
luminosity Lboi ~ 3.0/^0 and a circumstellar mass Mcnv~l-2 Mq (Shirley et al. 2000). 
Anglada et al. (1991) detected two radio continuum sources (VLAl and VLA 2) with 15" 
separation, both located within the error ellipse of the IRAS position. However, only VLA 2 
was found to be associated with dense gas in the cloud core (Girart et al. 1997). 

Figure 17a shows the SMA 1.3 mm dust continuum images of L723 VLA 2, in which two 
separated compact sources (labeled MMS 1 and MMS 2) are revealed. The angular separation 



39 



(3'/5 ± 0'/3) and mass ratio (1.0 ± 0.3) between the two sources derived in the SMA images are 
consistent with the resuhs derived earher in the OVRO 3 mm continuum observations (see 
Launhardt 2004). The two close compact sources are embedded in a common envelope seen 
in the SCUBA 850 /im observations. We note that the SMA data of L723 VIA 2 (including 
molecular lines results) have been published by Girart et al. (2009). Also note that VLA 
7 mm continuum observations discovered two components in source MMS 2 with an angular 
separation of ~ 1" (see Carrasco-Gonzalez et al. 2008), which would be unresolved by the 
SMA observations presented here. 



A. 22. B335 and L1157 

B335 is a well-studied, isolated Class protostar, and is widely accepted as one of the 
best protostellar collapse candidates (see, e.g., Zhou et al. 1993; Evans et al. 2005). It 
drives a large-scale bipolar outflow extending in the east-west direction (see, e.g., J0rgensen 
et al. 2007). L1157 is also a well-studied, isolated Class protostar, and is well-known 
for its prominent outflow. This outflow shows rich chemistry in various molecular lines 
observations (see, e.g., Bachiller et al. 2001; Arce et al. 2008), and is commonly considered 
the prototypical chemically active outflow. In the SMA 1.3 mm dust continuum observations, 
both B335 and L1157 remain as single sources (see Figure 17b & 17c), which is consistent 
with previous high angular resolution continuum observations toward the two sources (see, 
e.g., Harvey et al. 2003 for B335 and Chiang et al. 2010 for LI 157). 



A.23. L1251B 

L1251B is a small protostellar group located at 300±50pc (Kun & Prusti 1993). It is 
associated with an IRAS point source (IRAS 22376+7455), and was found to be associated 
with a CO outflow by Sato & Fukui (1989). Using OVRO and SMA (as well as single-dish 
telescopes), J. Lee et al. (2006; 2007) presented multiple-line observations toward L1251B, 
which shows very complex kinematics, including infall, rotation, and outflow motions. 

The SMA 1.3 mm dust continuum observations show four distinct sources embedded in 
the common envelope seen in the single-dish observations. Following J. Lee et al. (2007), 
the four sources are named L1251B IRSl, IRS2, SMA-N, and SMA-S (see Figure 17d). Both 
IRSl and IRS2 are associated with the infrared sources seen in the Spitzer images, and were 
classifled as Class O/I transition objects by J. Lee et al. (2006). Sources SMA-N and SMA-S 
have weak 1.3 mm dust continuum emission (see Figure 17d) and are not detected in the 
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Spitzer images, and are likely still in the early stages of the Class phase (see J. Lee et al. 
2007 for more details). 
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Fig. 1. — Distribution of angular resolutions and corresponding linear resolutions in this SMA 
survey. The dashed lines show the median values in the 1.3 mm continuum observations (3'.'1 
and 640 AU), the dotted lines show the median values in the 850 /xm continuum observations 
(I'.'O and 220 AU), while the solid lines show the median values for all observations (2'.'5 and 
600 AU). 
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Fig. 2. — Distribution of 1 a theoretical noises and 1 a measured noises in cleaned maps 
in this SMA survey. The dashed lines show the median values in the 1.3 mm continuum 
observations (1.0 and 2.3mJybeam~^), and the dotted lines show the median values in the 
850 /im continuum observations (1.7 and 5.5mJybeam^^). 
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Fig. 3. — (a) JCMT/SCUBA 850 /zm dust continuum image of L1448. Contour levels correspond to 1, 3, 5, 
8, 12, 15, and 20 a, then increase in steps of 15 cr (1 cr -40mJybeam-i). The SCUBA FWHM beam (- 14") 
is shown as a grey oval in the image, (b) SMA 1.3 mm dust continuum image of L1448N. Contour levels 
correspond to —3, 3, 6, 10, and 15cr, then increase in steps of lOcr (Icr ^ T.OmJybeam^^). (c) SMA 1.3mm 
dust continuum image of L1448 IRS2. Contour levels correspond to —3, 3, 6, 10, and 16 cr, then increase 
in steps of lOcr (Icr ^2.1 mJybeam^^). (d) SMA 1.3mm dust continuum image of L1448C. Contour levels 
correspond to —3, 3, 5, 9, and 15a, then increase in steps of 8(t (Icr ^4.0mJybeam^^). (e) SMA 850 /xm 
image of L1448C. Contour levels correspond to —3, 3, 6, 10, and 15 cr, then increase in steps of 10 cr (Icr 
^ 5.5inJy beam~^). The synthesized SMA beam is shown as a grey oval in the bottom right corner of each 
SMA dust continuum image. 
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Fig. 4. — SMA 1.3 mm dust continuum image of source PER 065, overlapped with the 
SCUBA 850 /im contours. The SMA contours (sohd hues) start from —3, 3 a, then increase 
in steps of 1 c (1 cr ~ 5.0 mJy beam"-*^). The synthesized SMA beam is shown as a grey oval 
in the bottom right corner. The SCUBA 850 /im contours (dotted lines) correspond to 75%, 
85%, and 95% of the peak emission (~0.3 Jybeam"^). 
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Fig. 5. — (a) JCMT/SCUBA 850 jum dust continuum image of the NCG 1333 central region. Contour levels 
correspond to 5, 10, 15, 20 and 30a, then increase in steps of 30 cr {la ^42mJybeam^^). The SCUBA 
FWHM beam is shown as a grey circle in the bottom right corner of the image, (b) SMA 1.3 mm dust 
contimumi image of SVS 13. Contour levels correspond to —3, 3, 6, 12, 20, and 30 cr, then increase in steps of 
20 cr (Icr ^2.5mJybeam~^). (c) Enlarged view of source SVS 13A. (d) SMA 1.3 mm dust continuum image 
of IRAS 2. Contour levels correspond to —3, 3, 6, 10, 15, 22, 30, and 40 cr, then increase in steps of 15 ct 
(Icr ~ 2.0mJybeam~^). (e) SMA 850 /xm dust continuum image of IRAS 2A. Contour levels correspond to 
—3, 3, 6, 10, 15, 22, and 30 cr, then increase in steps of lOcr (Icr ^4.5mJybcam^^). (f) SMA 1.3mm dust 
continuum image of IRAS 4. Contour levels correspond to —3, 3, 6, 12, 20, and 30 a, then increase in steps of 
15 cr (Icr ~ 12mJybeam~^). (g) SMA 850 ^m dust continuum image of IRAS 4A. Contour levels correspond 
to —3, 3, 6, 12, and 20 cr, then increase in steps of 10 cr (1 cr ^ 24mJybeam~^). The synthesized SMA beam 
is shown as a grey oval in the bottom right corner of each SMA dust continuum image. 
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Fig. 6. — (a) SMA 1.3 mm dust continuum image of IRAS 03282. The SMA contour levels 
correspond to —3, 3, 6, 10, 15, and 20 a, then increase in steps of 10 o" (1 cr ~ 1.8mJybeam~^). 
(b) SMA 850 /im dust continuum image of IRAS 03282. Contour levels correspond to —3, 
3, 5, 7, 11, 16, 22, and 30 a, then increase in steps of 10 a (la ~4.0mJybeam~^). The 
synthesized SMA beam is shown as a grey oval in the bottom right corner of the two SMA 
dust continuum images. 
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Fig. 7. — (a) JCMT/SCUBA 850 /xm dust continuum image of the Per-Bl central region, 
plotted on the Spitzer 24 /xm image. The SCUBA contour levels correspond to 10%, 20%, 
30%, 40%, 60%, 80%, and 90% of the peak emission (~ 2.5 Jybeam"!). The SCUBA FWHM 
beam is shown as a grey circle in the bottom right corner of the image, (b) SMA 1.3 mm 
dust continuum image of Per-Blc, overlaid with the SCUBA 850 ;um contours. The SMA 
contours (solid hues) correspond to —3, 3, 6, 10, 15, and 20 a, then increase in steps of 10 a 
(1 cr ~2.7mJybeam~^). The SCUBA 850 /im contours (dotted hues) correspond to 50%, 
70%, and 90% of the peak emission, (c) SMA 1.3 mm dust continuum image of Per-Blb, 
overlapped with the SCUBA 850 /im contours. The SMA contours (solid lines) correspond 
to —3, 3, 6, 10, 15, 20, 25, and 30 a, then increase in steps of 10 a {la ~4.2mJybeam^^). 
The SCUBA 850 /im contours (dotted lines) correspond to 30%, 50%, and 90% of the peak 
emission. The synthesized SMA beam is shown as a grey oval in the bottom right corner of 
the two SMA dust continuum images. 
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Fig. 8. — (a) SCUBA 850 /xm dust continuum image of the IC348 southern region. The 
SCUBA contour levels correspond to 3, 6, 10, 15, 20 and 25 a, then increase in steps of 
10 o" (1(7 ~40mJybeam~^). The SCUBA FWHM beam is shown as a grey circle in the 
bottom right corner of the image, (b) SMA 1.3 mm dust continuum image of IC348MMS, 
overlapped with the SCUBA 850 /im contours. The SMA contours (solid lines) correspond 
to —3, 3, 5, 7, 10, 15, 20, and 25 cr, then increase in steps of Gc {la ~4.0mJybeam^^). 
The SCUBA 850 /xm contours (dotted hues) correspond to 50%, 70%, and 90% of the peak 
emission of the IC 348 MMS core (~ 1.6 Jy beam-^). (c) SMA 850 /xm image of HH211 MMS. 
The SMA contours correspond to —3, 3, 6, 9, 12, 15, 20, and 25 a, then increase in steps 
of 5(7 (1(7 ~ 1.4mJybeam~^). The synthesized SMA beam is shown as a grey oval in the 
bottom right corner of the two SMA dust continuum images. 
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Fig. 9. — (a) SA-IA 1.3 mm dust continuum image of CB 17, overlapped with the SCUBA 850 /xm contours. 
The SMA contours (soUd hues) start at —3, 3cr, then increase in steps of Icr (1 cr~0.5mJybeam~^). The 
SCUBA contours (dotted hues) correspond to 50%, 70%, and 90% of the peak emission (~ 170mJybeam~^). 
(b) The same as Figure a, but for source IRAM 04191. The SMA contours start at —3, 3 a, then increase in 
steps of 1(7 (lcr--7mJybeam-i). The SCUBA 850 /zm contours correspond to 50%, 70%, and 90% of the 
peak emission (^0.37 Jybeam"^). (c) The same as Figure a, but for source L1521-IRS. The SMA contours 
correspond to —3, 3, and 5(7 (1 cr^ 1.5mJybeam^^). The SCUBA 850 /im contours correspond to 70%, 
85%, and 95% of the peak emission (^0.44 Jybcam^^). (d) The same as Figure a, but for source L1527. 
The SMA contours correspond to —3, 3, 6, 10, 15, 20, 30, and 40 a, then increase in steps of 15 u (Ict 
~ 1.5mJybeam~^). The SCUBA 850 y^m contours correspond to 50%, 70%, and 90% of the peak emission 
(^ 0.5 Jybeani^^). The synthesized SMA beam is shown as a grey oval in the bottom right corner of the 
each SMA image. 
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Fig. 10. — (a) SMA 1.3 mm dust continuum image of HH 114 MMS. The contour levels 
correspond to —3, 3, 6, 10, and 15 cr, then increase in steps of 10 o" (la ~5.0mJybeam~^). 
(b) SMA 880 /xm dust continuum image of HH114 MMS. The contour levels correspond to 
—3, 3, 6, 10, and 15 a, then increase in steps of 10 a (1 u ~ 16 mJy beam~^). The synthesized 
SMA beam is shown as a grey oval in the bottom right corner of the two SMA dust continuum 
images. 
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Fig. 11. — (a) SMA 850 /xm dust continuum image of 0MC3 MMS6, overlapped with the 
SCUBA 850 /xm contours. The SMA contours (sohd hues) correspond to —3, 3, 6, 10, 15, 20, 
and 25 a, then increase in steps of 8 a (1 cr ~ 20 mJy beam~^). The SCUBA 850 /im contours 
(dotted hues) correspond to 25%, 50%, 75%, and 95% of the peak emission (~27 Jybeam"^). 
(b) SMA 850 /xm dust continuum image of NGC2024FIR5. The SMA contours correspond 
to —3, 3, 6, 10, 15, and 20o", then increase in steps of lOo" (lo" ~ 14mJybeam~^). The 
synthesized SMA beam is shown as a grey oval in the bottom right corner of the two SMA 
dust continuum images. 
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Fig. 12. — (a) SMA 1.3 mm dust contimmm image of HH212 MMS (resolution 2'.'6 x 2'.'3). The SMA 
contours correspond to —3, 3, 5, 7, and lOcr, then increase in steps of 4cr [la ~'3.0mJybeam^^). (b) 
SMA 1.3mm dust continuum image (resolution C/.'S x (X'?); the contours correspond to —3, 3, 5, 7, and 
lOcr, then increase in steps of 5(7 (Icr ^ l.lmJybeam^^). (c) The SMA 850 /im dust continuum image 
(resolution 0'.'35 x 0'.'32); the contours correspond to —3, 3, 7, and 11 cr, then increase in steps of 4cr (Icr 
~ 2.8mJybeam^^). (d) The restored SMA 1.3 mm image with angular resolution of 0'.'32 x 0'.'32; the contours 
correspond to —3, 5, 8, 12, 16, 20, 25, and 30 ct, then increase in steps of 10 cr (Icr ~0.25mJybeam~^). The 
white crosses in Figures b & d show the position of the secondary continuum source derived from Figure c, 
while the stellar symbol in Figure c shows the position of the 1.4 mm dust continuum source detected by 
CodeUaet al. (2007). 
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Fig. 13. — (a) SMA 870 ;um dust continuum image of HH25 MMS, overlapped with the 
SCUBA 850 /im contours. The SMA contour levels (solid lines) correspond to —3, 3, 5, 8, 
and 12 a, then increase in steps of 8 ct (1 o" ~ 6.8 mJy beam"^). The SCUBA contours (dotted 
hues) represent the 50%, 75%, and 99% levels of the peak emission (~ 1.6 Jybeam"^). (b) 
SMA 1.3 mm dust continuum image of CG30, overlapped with the SCUBA 850 /xm contours. 
The SMA contours (solid lines) correspond to —3, 3, 6, 10, 15, 20, 25, and 30 a, then increase 
in steps of 10 a (la ~4.0mJybeam~^). The SCUBA contours (dotted lines) represent the 
50%, 75%, and 95% levels of the peak emission (~ 1.2 Jybeam"-*^). The synthesized SMA 
beam is shown as a grey oval in the bottom right corner of the two SMA dust continuum 
images. 
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Fig. 14. — (a) SMA 1.3 mm dust continuum image of B228. The contour levels correspond 
to —3, 3, 5, 7, and lOo", then increase in steps of Sc {\o ~ l.SmJybeam^-*^). (b) SMA 
1.1mm dust continuum image of B228. The contour levels correspond to —3, 3, 5, 7, 9, 11, 
14, and 18 o" (1 o" ~ 1.6 mJy beam~^). The synthesized SMA beam is shown as a grey oval in 
the bottom right corner of the two SMA images. 
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Fig. 15. — (a) SMA 1.3 mm dust continuum image of VLA 1623. The contour levels corre- 
spond to 5, 10, 15, 20, and 30o", then increase in steps of Scr (Ic ~ l.OmJybeam"^). (b) 
The enlarged view of the VLA 1623 East pair, (c) SMA 1.3 mm dust continuum image of 
Oph-MMS 126. The contour levels correspond to —3, 3, 4, and 5o" (lo" ~2.1mJybeam^^). 
(d) SMA 1.3 mm dust continuum image of IRAS 16293. The contour levels correspond to 
—3, 3, 5, 8, 12, 16, and 20o", then increase in steps of 5a {la ~ IQmJybeam^^). (e) The 
eSMA 850 ;um continuum image of source IRAS 16293 A. The contour levels correspond to 
—3, 5, and 10 o", then increase in steps of 5 0" {la ~20mJybeam~^). The synthesized SMA 
beam is shown as a grey oval in the bottom right corner of each SMA image. 
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Fig. 16. — SMA 850 ;um dust continuum image of L483. The contour levels correspond to 
—3, 3, 6, and 10 a, then increase in steps oi 5a {la ~3.2mJybeam~^). The synthesized 
SMA beam is shown as a grey oval in the bottom right corner. 
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Fig. 17. — (a) SMA 1.3 mm dust continuum image of L723, overlapped with the SCUBA 850 ^m contours. 
The SMA contours (soHd lines) correspond to —3, 3, 5, 8, 12, and 15 a levels and then increase in steps of 5 a 
(l(T--1.4mJybeam"i). The SCUBA 850 ^m contours (dotted lines) represent 30%, 50%, and 90% of the 
peak emission (^0.9 Jybeam^^). (b) Same as Figure a, but for source B335. The SMA contours correspond 
to —3, 3, 5, 8, 12, 16, and 20 cr levels and then increase in steps of 5cr (1 ct^ 2.9mJybeam~^). The SCUBA 
850 /im contours represent 30%, 50%, 70%, and 90% of the peak emission (^ 1.2 Jybcam^^). (c) Same as 
Figure b, but for source L1157, where la noise in the SMA image is ^4.8mJybeam^^ and the peak value 
in the SCUBA 850 /xm image is ~ 1.4 Jybeam"^. (d) SMA 1.3mm dust continuum image of L1251B. The 
SMA contours start at —3, 3(7 and then increase in steps of Icr (l(7~3.2mJybeam^^). 
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Fig. 18. — Top Left: Distribution of the 1.3 mm fluxes derived in this SMA survey. The me- 
dian value is 100 mJy. Top Right: Distribution of the 850 /xm fluxes (median value 430 mJy). 
Bottom Left: Distribution of the circumstellar gas masses derived from the 1.3 mm dust 
continuum observations (median value 0.06 M©). Bottom Right: Distribution of the cir- 
cumstellar gas masses derived from the 850 /im dust continuum observations (median value 
0.12 Mq). 
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Fig. 19. — Distribution of separations for protobinary systems in this SMA survey (see 
Table 6), binned in intervals of linear values of 1000 AU {top left) and 1600 AU {top right), 
and binned in intervals of log values of 0.25 {bottom left) and 0.40 {bottom right). The 
two dashed lines shown in the top figures mark the median linear resolution (600 AU) and 
median (projected) separation (1800 AU) in this survey, while the four dashed lines shown 
in the bottom figures mark the peak values in the distributions of separations for main 
sequence stars (30 AU; see Duquennoy & Mayor 1991) and pre-main sequence stars (60 AU; 
see Patience et al. 2002), and the median linear resolution and separation in this survey. We 
note that the 'separations' shown here are observed projected separations. 
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Fig. 20. — Observed multiplicity frequencies {left) and companion star fractions {right) 
with separations from 50 AU to 5000 AU for Class protostars (this work), Class I young 
stellar objects (data from Connelley et al. 2008a, b), and main sequence stars (data from 
Dquennoy & Mayor 1991). It must be noted that the SMA survey for Class protobinary 
systems is incomplete across the observed separations range, and the derived values should 
be considered as lower limits. 
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Fig. 21. — Histograms of circumstellar mass ratios for Class protobinary systems in this 
work (see Table 6), binned in intervals of 0.1 {left) and 0.2 {right). 



70 




- • • 

• - 

• • • - 

• • • 



4 6 8 10 12 14 16 18 

Angular Separation (arcsec) 



500 1000 1500 2000 2500 3000 3500 4000 4500 5000 

Projected linear Separation (AU) 



Fig. 22. — The circumstellar mass ratios plotted as a function of angular separations {left) 
and projected linear separations (right) for the Class protobinary systems in this work 
(see Table 6). The diagrams show that there is no obvious relation between mass ratios and 
separations for the protobinary systems in this sample. 
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Fig. 23. — Top: Temperature and density evolution at the center of a gravitationally collaps- 
ing cloud core obtained by Masunaga & Inutsuka (2000) in their radiation hydrodynamical 
calculation of spherically symmetric protostellar collapse. Each of these phases in the tem- 
perature evolution is characterized by a distinct value of the effective ratio of specific heats, 
7efj. The image is adopted from Andre et al. (2009). Bottom: A suggested picture of the 
sequential fragmentation for binary stars. From left to right: the IRAM-30m 1.3mm dust 
continuum image of Bok globule CB 246 (data from Launhardt et al. 2010), the SMA 1.3 mm 
dust continuum image of RCrA SMM lA (data from Chen & Arce 2010), and the eSMA 
850 /im image of IRAS 16293 A (this work), while the right blank panel shows a possible view 
with an angular resolution of 0'.'04 attained at the ALMA. 



